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Abstract 

Particle  image  velocimetry  (PIV)  has  been  carried  out  using  solid  carbon  dioxide 
(CO2)  particles  as  the  seed  material  to  continue  the  development  of  clean  seeding  for  use 
in  large-scale,  closed-circuit  tunnels.  Testing  occurred  in  two  wind  tunnels  at  subsonic 
and  supersonic  speeds  using  dry  ice  particles  generated  by  allowing  liquefied  C02  to 
expand  from  a  small  diameter  injector  tube  through  a  larger  diameter  shroud  tube. 
The  particles  were  injected  into  the  plenum  and  discrete  solid  particles,  suitable  for 
PIV  measurements,  were  present  in  the  test  section.  Data  on  particle  size  were  first 
collected  using  a  Malvern  particle  size  analyzer  for  three  sizes  of  injector  tubes,  two 
sizes  of  shroud  tubes,  and  two  different  types  of  shroud  tubes:  a  simple  tube  and  a 
static  mixing  tube.  The  injectors  using  the  static  mixing  shroud  tube  and  the  simple 
shroud  tube  were  each  used  in  the  adjustable  throat  supersonic  blowdown  wind  tunnel 
at  the  Air  Force  Institute  of  Technology  with  a  6  inch  by  6.5  inch  cross  section.  Particle 
size  results  for  these  two  configurations  suggested  that  the  static  mixing  shroud  tube 
decreased  the  Sauter  mean  particle  diameter  by  a  factor  of  three.  In  the  tunnel,  Mach  1.92 
flow  over  a  10°  ramp  was  produced  and  PIV  images  captured  particles  above  the  ramp 
both  upstream  and  downstream  of  the  oblique  shock  while  schlieren  imaging  provided 
insight  into  the  flow  conditions.  Both  the  velocities  far  upstream  and  far  downstream 
of  the  shock  closely  matched  expectations,  based  on  the  wind  tunnel  instrumentation. 
Particle  lag  for  the  flow  across  the  shock  was  quantified  for  the  two  cases,  and  despite  the 
substantial,  quantified  differences  in  particle  size  measured  at  the  shroud  tube  exit,  the 
results  for  both  shroud  tubes  were  generally  consistent  with  a  theoretical  response  of  a  2 
pm  particle.  Finally,  for  the  first  time  particles  were  injected  into  the  stilling  chamber  of 
the  Air  Force  Research  Laboratory’s  closed-circuit  Trisonic  Gas-dynamics  Facilty,  which 
has  a  24  inch  by  24  inch  cross-section  at  three  subsonic  speeds  and  four  stagnation 
pressures.  PIV  was  successfully  carried  out  in  each  case.  Measured  streamwise  velocities 
matched  expected  velocities  within  a  few  percent  based  on  tunnel  instrumentation,  and 


IV 


freestream  turbulence  was  found  to  be  less  than  2%  in  most  cases.  These  results  suggest 
that  PIV  using  CO2  particles  may  be  robustly  implemented  in  this  closed-circuit  wind 
tunnel  without  risk  of  contaminating  the  tunnel. 
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I.  Introduction 

AVIATION  has  made  quick  and  efficient  world  travel  possible,  has  given  a  new  di¬ 
mension  to  the  art  of  war,  and  has  become  the  stepping  stone  into  the  space 
age.  Improving  this  industry  has  motivated  the  ever  increasing  study  of  fluid  mechanics, 
allowing  the  advancement  of  aircraft  and  rockets  alike.  It  stands  to  reason  that,  with 
such  a  vast  impact  on  these  industries,  improving  the  methods  and  efficiency  with  which 
this  understanding  is  achieved  also  has  its  merits.  In  modern  history,  theory,  observation 
and  experimentation  have  been  the  predominant  means  to  make  advancements  in  the 
field.  Though  there  is  limited  evidence  of  his  experiments  in  flight  actually  being  carried 
out,  the  observations,  drawings  and  notes  of  Leonardo  da  Vinci  showed  an  analytical 
process  to  develop  a  way  for  humans  to  mimic  the  flight  of  birds  [1].  Since  his  time, 
countless  individuals  have  tackled  issues  dealing  with  the  mechanics  of  fluids,  including 
Newton,  Bernoulli,  Euler,  Navier,  Stokes,  Reynolds,  and  Prandtl,  just  to  name  a  few. 

The  fundamental  theories,  concepts  and  equations  developed  by  these  and  many 
others  are  the  basis  from  which  human  flight  was  made  possible,  as  proven  by  the  Wright 
brothers  in  the  early  1900s.  The  Wrights,  however,  experienced  setbacks  with  lower  than 
expected  performance  of  their  original  gliders  that  could  have  halted  their  progress  if 
they  were  without  ingenuity  and  perseverance  [2],  In  late  1901  they  designed  and  built 
a  wind  tunnel  to  test  different  wing  designs,  ultimately  providing  the  data  necessary 
to  move  forward  with  their  aircraft.  As  it  was  then,  today  it  is  still  necessary  to  test 
items  of  interest  in  facilities  such  as  wind  tunnels  because,  even  with  the  advent  and 
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rise  of  computer  modeling,  there  are  still  many  problems  in  fluid  motion  that  require 
experimental  data. 

1 . 1  Motivation 

Many  experimental  methods  exist  today  that  can  be  used  to  gain  insight  into 
an  area  of  interest,  each  with  its  advantages  and  disadvantages.  One  of  these  is  particle 
image  velocimetry  (PIV)  which  is  based  on  the  simple  relation  that  velocity  is  the  distance 
traveled  per  unit  of  time.  If  particles  can  be  placed  in  a  fluid  under  the  right  conditions, 
the  particles  will  move  at  the  same  speed  as  the  fluid.  If  two  sequential  images  of  a  particle 
are  captured  with  a  known  time  increment  between  the  two,  the  distance  traveled  by  the 
particle  can  be  deduced  by  comparing  the  first  image  to  the  second.  Coupled  with  the 
knowledge  of  the  time  difference  between  the  two  images,  the  speed  can  be  determined. 
Raffel,  et.  al.  make  the  analogy  between  this  and  children  throwing  pieces  of  wood  in 
a  stream  to  observe  them  float  downstream  [3].  A  practical  application  of  this  simple 
concept  has  been  applied  by  sailors  for  centuries  when  a  floating  object  such  as  a  log  or 
buoy  tied  to  a  rope  would  be  tossed  overboard.  Equally  spaced  knots  would  be  tied  in 
the  rope  and  as  the  rope  unraveled  overboard,  the  number  of  knots  would  be  counted 
and  timed,  likely  with  a  sand  glass.  When  the  sand  ran  out  the  number  of  knots  would 
be  tallied  and  the  ship’s  speed  was  determined.  Hence  the  nautical  term  “knot”  was 
coined  to  describe  speed  [4], 

As  simple  as  this  concept  is,  the  application  to  aerodynamics  and  other  fields  of 
study  can  be  quite  complex,  as  discussed  in  chapter  2,  due  in  part  to  the  presence  of  mul¬ 
tiple  particles  within  the  held  of  view.  Once  these  complexities  are  understood,  however, 
a  great  deal  of  information  can  be  obtained  about  the  fluid  how  of  interest.  With  planar 
velocity  data,  one  can  determine  boundary  layer  effects,  turbulence  characteristics,  eddy 
formations,  momentum  thickness,  and  so  on.  The  applications  of  PIV  in  aerodynamics 
and  water  hows  continue  to  grow  [3].  Consequently,  there  is  interest  in  refining  this 
technique  to  improve  its  accuracy  as  well  as  its  practicality. 

One  of  the  problems  presented  by  PIV  is  that  introducing  particles  in  tunnels  can 
be  inconvenient  and  costly.  One  facility,  the  Air  Force  Research  Laboratory’s  (AFRL) 
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Trisonic  Gas-dynamics  Facility  (TGF),  requires  additional  maintenance  to  clear  the  tun¬ 
nel  of  the  particles  used  for  PIV  measurements.  The  TGF  is  a  closed-circuit  supersonic 
wind  tunnel,  meaning  that  particulates  introduced  into  the  flow  for  PIV  persist  either 
in  the  flow  or  accumulate  on  surfaces  within  the  tunnel.  In  one  recent  case,  the  use  of 
an  oil-based  smoke  to  seed  the  air  flow  caused  a  residue  to  build  up  on  the  tunnel  walls 
and  test  section  glass,  necessitating  substantial  cleaning  after  the  test.  Also,  residue 
from  oil-based  or  solid  seed  material,  such  as  titanium  dioxide,  can  spoil  measurements 
executed  using  pressure-sensitive  paint.  Additional  maintenance  means  additional  cost 
and  increased  downtime  of  the  tunnel.  With  this  type  of  impact,  the  practicality  of 
using  PIV  in  the  TGF  and  other  similar  facilities  is  reduced  until  there  exists  a  method 
to  cleanly  seed  the  flow.  Many  examples  of  oil-  and  water-based  aerosolized  seeding  are 
available,  including  recent  studies  at  the  Arnold  Engineering  Development  Center’s  16T 
wind  tunnel,  showing  that  if  these  seed  materials  can  be  replaced  with  clean  seeding 
methods,  many  facilities  could  benefit  [5,  6].  To  this  end,  the  past  several  years  have 
seen  the  introduction  and  advancement  of  using  solid  carbon  dioxide  (C02),  or  dry  ice, 
to  seed  the  flow.  C02  can  exist  as  a  solid  and  then  sublimate,  leaving  no  trace  or  residue 
other  than  the  naturally  occurring  inert  gas  that  mixes  into  the  surrounding  air.  Thus, 
if  particles  of  dry  ice  can  be  inserted  into  the  air  flow,  they  can  persist  long  enough 
to  pass  through  the  test  section  and  then  sublimate,  negating  any  need  for  additional 
maintenance  on  the  wind  tunnel  due  to  PIV  testing  alone. 

1.2  Research  Focus 

The  selection  of  particle  type  is  critical  for  successful  PIV  measurements,  not  only 
because  of  the  impact  to  the  facilities  using  PIV,  but  for  the  measurements  themselves. 
It  is  necessary  to  select  particles  with  the  proper  characteristics  that  make  for  accurate 
measurements,  including  size,  density,  quantity,  and  light  scattering  characteristics.  The 
research  presented  here  focuses  on  correctly  sizing  C02  particles  and  determining  meth¬ 
ods  to  increase  the  number  of  particles  to  provide  sufficient  data  for  accurate  velocity 
vector  mapping.  Particles  must  be  small  enough  to  accurately  track  the  fluid  in  ques¬ 
tion  while  still  being  large  enough  to  scatter  ample  light  to  be  captured  by  a  camera 
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sensor.  The  particle  density  should  be  as  close  to  the  fluid  density  as  possible  so  that 
the  effect  clue  to  gravity  and  the  response  time  clue  to  fluid  accelerations  are  minimized. 
Lastly,  there  needs  to  be  sufficient  particle  coverage  over  the  field  of  view  to  maximize 
the  resolution  of  the  velocity  vector  field. 

This  research  also  attempts  to  scale  the  seeding  method  by  injecting  particles  into 
the  TGF  and,  for  the  first  time,  collecting  PIV  data  using  CO2  seeding.  Previous  work 
by  Greene,  DeLapp,  McNiel,  and  Peltier,  under  the  direction  of  Dr.  Mark  Reeder,  honed 
the  methods  for  C02  particle  generation  and  its  viability  for  PIV  application  [7,  8,  9]; 
this  research  refines  their  methods  with  the  following  goals: 

•  Increase  particle  count  to  improve  spatial  resolution  of  measurements  in  the  AFIT 
supersonic  wind  tunnel 

•  Assess  methods  for  controlling  particle  size 

•  Characterize  particle  response  to  an  oblique  shock  wave  using  different  injectors 

•  Demonstrate  PIV  in  the  TGF 

Previously,  the  injectors  used  for  seeding  were  varied  to  assess  the  size  distribution  on 
the  resulting  particles.  Factors  such  as  the  injection  orifice  size,  injector  shroud  size,  and 
injection  point  were  all  found  to  be  important.  Several  injection  methods  and  injection 
points  were  tested  to  find  a  suitable  solution  to  the  seeding  problem.  However,  attempts 
to  scale  the  injection  method  for  use  in  a  larger  wind  tunnel  had  limited  success  in 
achieving  sufficient  particle  coverage  [10,  11].  Thus  the  research  presented  in  the  following 
chapters  discusses  the  attempt  to  improve  upon  this  along  with  the  goals  listed  above. 
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II.  Background 


IN  this  chapter,  background  information  is  provided  on  the  systems  used  for  this  re¬ 
search.  First,  the  basics  of  PIV  are  covered,  discussing  the  components  that  are 
necessary  to  accomplish  PIV.  Emphasis  is  placed  on  flow  seeding  because  of  its  impor¬ 
tance  in  this  research.  Recent  research  in  the  held  of  clean-seeding  is  discussed.  Next, 
the  wind  tunnel  and  the  mathematical  relations  that  will  be  relevant  for  calculations  in 
this  research  are  discussed.  Finally,  the  basic  schlieren  setup  is  discussed  as  it  pertains 
to  this  research. 

2.1  PIV  Overview 

Advancements  in  technology  over  the  past  decades  have  allowed  PIV  to  become  a 
viable  and  common  measurement  technique.  Correlating  particle  motion  between  differ¬ 
ent  images  is  a  hurdle  that  has  been  relatively  easy  to  overcome  with  the  introduction  of 
digital  charge  coupled  device  (CCD)  camera  technology  and  computing  power  enhance¬ 
ments.  Figure  2.1  displays  the  basic  components  of  PIV,  namely  the  seeding  material, 
light  source,  imaging  optics  and  imaging  sensor.  Each  of  these  components  is  briefly 
described  below  but  the  interested  reader  is  referred  to  Raffel,  et.  al.,  for  an  in  depth 
review  of  PIV  [3], 

2.1.1  Light  Source.  Light  is  necessary  to  illuminate  the  particles  in  the  flow 
sufficiently  to  be  detected  by  the  camera  sensor.  For  planar  PIV  as  shown  in  Figure  2.1, 
lasers  are  typically  used  because  they  emit  monochromatic  light  with  high  energy  density 
that  is  easily  manipulated  through  a  careful  selection  of  optics  to  produce  a  light  sheet 
that  illuminates  the  full  held  of  view  of  the  camera  [3].  A  common  type  of  laser,  also 
used  for  this  research,  is  the  neodymium-doped  yttrium-aluminum-garnet  (ND:YAG) 
laser  which  uses  two  resonators  that  produce  light  at  a  wavelength  of  A  =  1064  nm 
that  is  then  frequency  doubled  to  get  a  wavelength  of  A  =  532  nm  in  the  visible  green 
spectrum.  Solid  state  lasers  such  as  this  can  produce  high  repetition,  low  duration  pulses 
of  energy  making  them  ideal  for  PIV  applications. 


5 


Mirror 


Figure  2.1:  Principle  components  of  PIV. 

[3] 

Selection  of  a  light  source  also  depends  on  the  scattering  characteristics  of  the 
particles  used  for  seeding.  If  the  particles  used  are  smaller  than  the  wavelength  of  light, 
the  light  is  scattered  according  to  Rayleigh’s  theory  where  the  view  angle  is  less  important 
because  the  light  is  scattered  more  evenly  in  all  directions.  The  energy  incident  on  the 
sensor  is  accordingly  less  so  detection  of  the  particles  becomes  more  difficult.  For  particles 
larger  than  the  wavelength  of  the  light,  the  light  is  scattered  according  to  Mie  theory 
where  there  is  a  preponderance  of  light  scattered  in  the  same  direction  as  the  incumbent 
light  energy  and  less  scattered  backwards  and  to  the  sides.  Figure  2.2  compares  light 
scattering  for  Rayleigh  and  Mie  theory.  Table  2.1  lists  most  of  the  typical  seed  materials 
for  air  flow  and  their  sizes.  It  is  apparent  that  for  PIV  applications  using  Nd:YAG 
lasers,  Mie  scattering  applies  because  the  particles  are  typically  larger  in  diameter  than 
the  wavelength  of  the  light. 

2.1.2  Flow  Seeding.  Table  2.1  shows  common  seeding  types  used  for  PIV 
measurements.  The  particles  chosen  for  a  specific  test  must  be  able  to  scatter  sufficient 
light,  but  must  also  respond  quickly  to  the  dynamics  of  the  flow  being  studied.  Larger 
particles  will  scatter  more  light  than  smaller  ones  of  the  same  material,  but  a  larger 
particle  has  more  momentum  at  the  same  velocity  as  the  smaller  one,  meaning  that 
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Rayleigh  scattering 


Mie  scattering 


Figure  2.2: 


Rayleigh  and  Mie  scattering  characteristics. 


[12] 


Table  2.1:  Common  seeding  materials  for  gas  flows. 


Type 

Material 

Mean  diameter  in  pm 

Solid 

Polystyrene 

0.5 

10 

Alumina  AI2O3 

0.2 

-  5 

Titania  TiOa 

0.1 

-  5 

Glass  micro-spheres 

0.2 

-  3 

Glass  micro-balloons 

30 

-  100 

Granules  for  synthetic  coatings 

10 

-  50 

Dioctylphathalate 

1 

-  10 

Smoke 

<  1 

Liquid 

Different  oils 

0.5 

-  10 

Di-ethyl-hexyl-sebacate  (DEHS) 

0.5 

-  1.5 

Helium-filled  soap  bubbles 

1000 

-  3000 
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more  force  will  be  required  to  affect  the  particle’s  momentum.  For  a  spherical  particle 
of  constant  density,  the  mass  increases  with  the  cube  of  the  radius,  r3,  whereas  the 
viscous  forces  only  act  on  the  surface  area  of  the  sphere  which  increases  with  r2.  Thus, 
the  momentum  increase  with  particle  size  is  greater  than  the  available  aerodynamic 
force  acting  on  the  particle  to  change  the  momentum,  thus  larger  particles  will  lag  the 
dynamics  of  the  air  flow  more  than  the  smaller  ones.  Various  forms  of  the  Bassett- 
Boussinesq-Oseen  (BBO)  equation  of  motion  for  a  spherical  particle  in  a  viscous  fluid 
have  been  noted  in  the  literature  to  characterize  particle  response  [9].  Melling  gives  this 
equation  as  [13]: 


vr  d3  d  Up 
-PP- 
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If  the  density  of  the  particle  is  much  larger  than  the  density  of  the  fluid,  as  is  the  case 
with  many  seeding  methods  in  air,  the  BBO  equation  can  be  simplified  significantly  to 
equation  2.2  [13].  The  CO2  density  is  assumed  to  be  1.18  g/cm3  based  on  previous  work, 
so  that  this  simplification  is  valid  [10]. 


dvp 

dt 


-C(vp  -  U ) 


(2.2) 


In  equation  2.2  vp  is  the  particle  velocity,  U  is  the  fluid  velocity,  and  C  is  a  characteristic 
frequency  of  the  particle  motion  in  terms  of  the  drag  coefficient  and  Reynolds  number. 
For  Stokes’  resistance,  the  drag  coefficient  is  given  by  Co  =  24 / Rep,  yielding  the  following 
result  for  C: 


18/i 

Pvdl 


(2.3) 


For  very  small  particles  or  rarefied  flows,  a  Knudsen  number  correction  can  be  applied 
to  the  drag  coefficient  to  account  for  the  ratio  increase  between  the  mean  free  path  of 
molecules  in  the  fluid  to  the  particle’s  diameter.  The  drag  coefficient  would  then  be 
given  by  Co  =  24  [Rep  (1  +  Knp)]~l  [13].  In  Greene’s  research,  he  compared  his  results 
using  a  drag  coefficient  with  and  without  the  Knudsen  number  correction  and  found  that 


the  Knudsen  number  had  very  little  effect  [10].  For  this  reason,  the  Knudsen  number 
correction  is  not  used  in  the  current  research. 

The  Stokes’  number  is  often  defined  as 

PpdpUCgiip 

=  18  fiL 

where  U  is  the  velocity,  Csnp  is  a  constant  used  to  account  for  velocity  slip  in  rarefied  flow 
fields,  and  L  is  a  characteristic  length  scale  of  the  flow.  Here,  Csup  is  approximately  unity. 
In  a  turbulent  flow,  L  is  generally  taken  to  be  either  the  integral  scale  or  Kolmogorov 
scale.  In  this  study,  however,  the  shock  thickness  is  used  to  characterize  the  flow.  Since 
this  thickness  is  very  small,  the  Stokes’  number  would  increase  toward  infinity.  So  rather 
than  utilizing  Stokes’  number  to  characterize  the  particle  response,  one  goal  of  this 
project  is  to  characterize  the  particles  by  another  means. 

2.1.3  Image  Acquisition  and  Correlation.  Image  acquisition  in  PIV  is  accom¬ 
plished  primarily  with  digital  cameras.  CCD  and  complementary  metal-oxide  semicon¬ 
ductor  (CMOS)  sensors  convert  light  energy  into  electrical  energy  that  is  recorded  and 
processed  to  produce  images.  While  CCD  cameras  typically  have  a  higher  resolution  than 
CMOS,  CMOS  sensors  are  more  sensitive  to  light,  making  them  better  suited  for  low 
light  conditions.  In  either  case,  the  shutter  is  typically  open  longer  than  the  pulse  of  laser 
energy  so  low  ambient  light  conditions  are  favorable  for  PIV.  This  reduces  background 
noise  and  reflections  in  the  images  acquired.  In  the  past,  PIV  images  were  collected 
on  film,  sometimes  using  a  single  exposure  to  capture  both  laser  pulses.  Today,  digital 
cameras  are  able  to  capture  two  images  in  rapid  succession,  allowing  the  first  image  to 
capture  one  laser  pulse  and  the  second  image  to  capture  the  next  laser  pulse.  Digital 
processing  allows  for  highly  automated  correlations  of  the  two  images  using  prepackaged 
software  programs.  A  cross-correlation  method  splits  the  image  into  small  subregions, 
called  interrogation  regions  (IR),  and  then  uses  a  fast  Fourier  transform  (FFT)  on  each 
IR.  The  signal  strength  for  a  single  IR  is  compared  between  the  first  and  second  image 
to  determine  a  mean  displacement  between  the  two  peaks.  The  correlation  produces  a 
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vector  for  that  IR,  and  then  the  next  IR  is  processed.  This  process  is  shown  in  Figure 
2.3  below. 


System 

input 


Image  subsampling 
at  position  (i  j) 


Fast  Fourier 
transformation 


Correlation  in  the 
spatial  frequency 
domain 


Inverse  Fourier 
transformation 


Peak  detection  and 
subpntd  interpolation 


Conversion 
to  velocity 


System  output 


Image  I:  Image  2: 
t-t0  l=to+At 


Window  function 
(optional) 


Filter  function 
(optional) 


Data  file 


Figure  2.3: 


Diagram  of  the  cross  correlation  process. 
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2.2  Discrete  COo  Particle  Formation  and  Characterization 

The  process  of  forming  discrete  CO2  particles  has  been  well  documented  in  work 
preceding  this  research  [7,  8,  9].  Pressurized  liquid  CO2  flows  through  an  expansion 
nozzle  at  which  point  some  of  the  liquid  evaporates  and  cools  rapidly.  The  cool  gas 
drops  the  temperature  of  the  remaining  liquid  spray,  freezing  it  into  solid  particles  that 
can  agglomerate  to  form  variously-sized  discrete  particles.  Initially,  this  method  for 
forming  particles  was  applied  to  the  cleaning  of  substrates,  as  patented  by  Whitlock, 
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Weltmer  and  Clark  [15].  In  2006  DeLapp,  Reeder,  Grafton  and  Goss  proposed  using 
these  particles  to  seed  air  flow  for  PIV  measurements  [7].  Their  work  initially  used  a 
commercial  off  the  shelf  dry  ice  cleaning  system  called  the  Sno-Gun  II  to  generate  CO2 
particles.  The  particle  stream  was  passed  through  a  Malvern  particle  size  analyzer  using 
different  nozzle  sizes,  resulting  in  sizes  typically  ranging  from  5  to  15  pm. 

McNiel  continued  and  built  upon  the  work  of  DeLapp  and  assessed  whether  or 
not  discrete  particles  could  be  generated  and  inserted  into  the  air  flow  of  a  wind  tunnel 
for  PIV  measurements  [8].  McNiel  studied  two  different  injection  methods:  a  capped, 
multi-port  shroud  tube  and  a  simple,  uncapped  shroud  tube.  He  also  used  these  two 
injection  methods  in  two  different  locations:  the  stilling  chamber  and  the  sidewall  of  the 
converging-diverging  nozzle.  McNiel  found  that  using  a  multi-port  injection  method  was 
unsuccessful  in  forming  discrete  particles  from  either  location  and  instead,  gaseous  C02 
clouds  were  present  in  the  test  section.  Using  the  simple  shroud  tube,  he  successfully 
generated  discrete  particles  that  persisted  through  the  test  section  at  the  nozzle  injection 
point  [8]. 

The  work  of  Peltier,  in  addition  to  McNicl’s  work,  explored  many  combinations 
of  injection  location  and  method,  moving  away  from  the  Sno-Gun  apparatus  and  using 
a  combination  of  small  and  large  tubes  to  constrict  and  expand  the  liquid  C02  in  a 
controlled  fashion  [8].  Greene  used  the  results  of  this  previous  work  and  studied  the 
effects  of  varying  the  size  of  the  feed  tube  and  the  shroud  tube  diameters,  as  well  as  the 
addition  of  a  coflow  of  air  to  the  C02  spray  [9] .  Another  variable  Greene  looked  at  is  the 
effect  insulating  the  tubes  had  on  particle  size.  His  work  using  the  particle  size  analyzer 
found  that  lengthening  the  shroud  tube,  increasing  the  shroud  tube  inner  diameter  (ID), 
decreasing  the  feed  tube  ID,  and  insulating  the  shroud  tube  each  effectively  increased  the 
size  of  particles  generated  [10].  Conversely,  particle  size  would  decrease  if  any  of  these 
changes  were  reversed.  The  addition  of  a  coflow  of  air  decreased  the  particle  size  but 
the  net  effect  was  a  reduction  in  the  particle  number  density  in  the  test  section.  Greene 
focused  on  injecting  particles  and  collecting  PIV  data  in  the  AFIT  2.5  inch  by  2.5  inch 
pressure-vacuum  supersonic  wind  tunnel  with  great  success.  His  attempt  to  scale  the 
injection  system  up  to  the  larger  6  inch  by  6.5  inch  supersonic  blowdown  tunnel  was  met 
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with  some  success.  However,  particle  number  density  was  lower  than  that  observed  in 
the  smaller  tunnel.  He  attributed  this  to  the  increase  in  air  mass  flow  rate  for  the  larger 
tunnel  and  the  increased  distance  between  the  injection  point  and  the  test  section,  which 
in  turn,  increased  the  required  residence  time  for  the  particles  to  persist  through  the  test 
section  [9]. 

2.3  Particle  Size  Analysis 

The  behavior  of  light  when  it  interacts  with  particles  allows  a  passive  means  for 
determining  particle  size.  Light  impinging  on  a  particle  will  be  scattered  in  a  well  known 
manner  depending  on  the  wavelength  of  the  light  and  the  size  of  the  particle.  If  the 
wavelength  is  known  and  an  array  of  sensors  are  used  to  detect  the  scattered  energy 
from  a  particle,  then  the  size  of  the  particle  can  be  deduced.  This  method  is  exploited 
by  particle  size  analyzers  which  shine  a  laser  across  a  known  distance  to  a  lens.  The  lens 
focuses  the  light  on  the  center  of  a  ring  detector.  When  something  such  as  a  particle 
disrupts  the  path  of  the  light,  the  light  will  diffract  and  impact  a  different  detector.  To 
arrive  at  a  particle  size,  however,  the  shape  is  assumed  spherical  because  a  sphere  is 
the  only  shape  that  can  be  uniquely  defined  by  a  single  number,  namely  its  diameter. 
For  particles  that  are  not  spherical,  the  true  size  will  be  different  than  that  found  by 
size  analyzers,  which  instead  provide  an  equivalent  spherical  size.  Greene  documented 
that  the  Sauter  mean  diameter,  D32,  was  a  useful  product  in  arriving  at  a  mean  particle 
diameter  [10].  It  is  calculated  by  summing  the  cube  of  the  diameters  and  then  dividing 
this  by  the  sum  of  the  diameters  squared  to  obtain  a  mean  particle  diameter,  as  seen  in 
equation  2.5  [16]. 

N 

D32  =  ^ -  (2.5) 

XLi-0? 

i= 1 

D32  is  the  mean  diameter  of  an  equivalent  sphere  that  has  the  same  volume/surface  area 
ratio  as  the  particles  being  measured.  In  general,  the  Sauter-mean  diameter  is  commonly 
used  in  experimental  studies  [17]. 
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2-4  High  Speed  Wind  Tunnel  Testing 

Wind  tunnels  are  used  to  simulate  actual  flight  in  a  controlled  environment,  thereby 
providing  a  means  for  acquiring  data  without  subjecting  the  test  to  undue  risks.  Uni¬ 
form  air  flow  at  subsonic,  transonic,  and  supersonic  speeds  can  be  achieved  at  variable 
conditions  for  easy  prediction  of  aerodynamic  characteristics  in  a  real  environment. 


2-4-1  Tunnel  Design.  All  tunnels  operate  by  creating  a  pressure  gradient 
across  a  control  volume  to  induce  air  flow  from  high  to  low  pressure.  High  pressure  can 
be  created  by  using  fans,  compressors,  or  a  supply  of  compressed  air  that  is  released 
into  the  lower  ambient  pressure  through  the  tunnel.  Some  systems  combine  this  with  a 
vacuum  system  to  further  increase  the  pressure  differential.  Systems  that  use  supplies  of 
compressed  air  and/or  vacuums  to  achieve  the  flow  cannot  run  continuously  because  they 
exhaust  their  supply  of  compressed  air,  whereas  systems  that  use  compressors  can  run 
continuously.  To  achieve  supersonic  flow,  it  is  necessary  to  design  a  compressor-diffuser 
nozzle  through  which  the  air  must  flow.  When  the  pressure  differential  is  large  enough, 
the  flow  will  choke  at  the  throat  of  the  nozzle  causing  a  normal  shock  to  form.  The  air 
then  expands  and  accelerates  based  on  the  ratio  of  the  area  of  the  test  section  to  the  area 
of  the  throat.  The  volume  upstream  of  the  throat  is  often  referred  to  as  the  stilling  or 
stagnation  chamber,  and  if  the  air  properties  are  known  inside  this  region,  the  properties 
of  the  flow  through  the  tunnel  can  easily  be  calculated  using  the  following  relations: 


P 

Po 

P_ 

Po 

T 

To 


(2.6) 


(2.7) 


(2.8) 
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2-4-2  Mass  Flow.  The  conservation  of  mass  can  be  used  to  determine  the  mass 


flow  of  air,  rh,  through  the  wind  tunnel  using  equation  2.9. 


m  =  pU  A 


(2.9) 


The  density,  p,  can  be  found  using  equation  2.7,  and  the  area,  A,  is  known  from  the 
tunnel  geometry  so  all  that  is  needed  is  the  velocity,  U .  Using  pressure  transducers  in  a 
wind  tunnel,  the  pressure  ratio  Pq/P  can  be  measured  so  that  equation  2.6  can  be  solved 
for  Mach  number: 


M 


(2.10) 


Since  the  Mach  number  is  defined  as  M  =  U /a  with  the  speed  of  sound  defined  as 
a  =  y/'yRT,  the  velocity  can  be  calculated.  Thus  the  mass  flow  rate  of  air  can  be  found 
using  pressure  and  temperature  measurements  in  the  tunnel. 


2-4-3  Isentropic  Relations.  Previous  work  by  Greene  looked  at  a  10°  ramp  in 
the  AF1T  2.5  inch  by  2.5  inch  pressure-vacuum  wind  tunnel,  so  for  similarity,  the  same 
will  be  done  for  this  research  in  the  larger  6  inch  by  6.5  inch  blowdown  wind  tunnel. 
The  ramp  with  known  deflection  angle,  9 ,  will  cause  an  oblique  shock  wave  at  an  angle 
(3.  If  the  speed  of  the  flow  is  known,  then  the  well  known  6  —  /3  —  M ach  relation  given 
by  equation  2.11  can  be  used  to  determine  the  shock  angle. 


tan  9 


2  cot  f3 


M 2  sin2  (3  -  1 
M2  (7  +  cos  2/3)  +  2 


(2.11) 


Knowledge  of  both  the  upstream  Mach  number  and  the  shock  angle  provides  the  nec¬ 
essary  information  to  determine  the  velocity  normal  to  the  shock  upstream  and  down¬ 
stream  of  the  shock  wave  using  isentropic  relations.  First,  the  Mach  number  normal  to 
the  shock  is  simply  computed  and  then  used  to  determine  the  Mach  number  downstream 
and  normal  to  the  shock. 

Mnl  =  M 1  sin  [3  (2-12) 
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(2.13) 


Mi,  = 


i  +  (¥Xi 


"•2  7 Ml  -  ^ 

Finally  the  true  Mach  number  downstream  of  the  shock  can  be  calculated  using  the 


following: 


M2  = 


M, 


n,  2 


sin  {/3  -  9) 


(2.14) 


2.5  Schlieren  Imaging 

“Light  propagates  uniformly  through  homogeneous  media”  but  small  disturbances 
due  to  temperature  variations,  turbulence,  convection,  weather  phenomena  or  shock 
waves  in  air  affect  the  density  of  the  air,  thereby  changing  its  refractive  index  [18]. 
Air  has  a  refractive  index  of  approximately  n  =  1.0003  but  localized  disturbances  can 
change  this  slightly,  varying  with  density  linearly  according  to  the  simple  relationship 
from  Settles  [18]: 

n  —  1  =  kp  (2.15) 

The  Gladstone-Dale  coefficient,  k,  is  about  0.23  cm3/g  for  air  at  standard  conditions 
[18].  Since  the  speed  of  light  through  a  media  other  than  a  vacuum  is  proportional  to 
the  index  of  refraction  according  to  equation  2.16,  then  a  change  of  the  refractive  index 
due  to  a  density  gradient  affects  the  speed  light  traverses  through  that  media,  causing  a 
phase  shift. 

A  v=-  (2.16) 

n 

Here,  A  is  the  wavelength,  v  is  the  frequency,  and  c  is  the  speed  of  light.  An  example 
of  this  is  easily  observed  when  looking  at  an  object  through  the  heat  rising  off  a  flame. 
These  gradient  disturbances  are  known  as  schlieren  and  they  bend  the  light  slightly, 
typically  too  little  to  be  seen  by  the  naked  eye.  If,  however,  collimated  light  is  passed 
through  the  schlieren  and  then  focused  to  a  point  by  a  lens,  the  schlieren  will  change 
the  path  of  the  light  through  the  lens,  causing  a  change  in  the  point  at  which  the  light 
focuses.  By  placing  a  knife  edge  at  the  true  focal  point  of  the  lens,  portions  of  the  bent 
light  are  blocked  by  the  knife  edge,  causing  a  shadow  to  form  against  an  otherwise  bright 
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background.  This  is  demonstrated  in  the  simple  diagram  from  Settles,  seen  in  Figure 
2.4. 


lens  lens  screen 


Figure  2.4: 


Diagram  of  a  simple  schlieren  system  with  a  point  light  source. 


[18] 


A  popular  arrangement  is  the  z-type  Herschellian  system  shown  in  Figure  2.5  be¬ 
cause  mirrors  can  provide  a  larger  field  of  view  than  lenses  of  equal  cost  [18].  Light  is 
condensed  through  a  slit  positioned  at  the  focal  point  of  a  parabolic  mirror  such  that 
the  reflected  light  is  collimated.  The  light  then  passes  through  the  area  of  interest  and 
is  bounced  off  another  mirror  and  is  focused  onto  the  knife  edge,  blocking  part  of  the 
light.  A  camera  captures  the  light  which  passes  through  as  well  as  the  shadows  caused 
by  any  disturbances  in  the  path  of  the  light. 


Figure  2.5: 


Diagram  of  a  Herschellian  schlieren  system. 


[18] 
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III.  Methodology 


This  chapter  discusses  the  test  equipment  and  test  procedures  used  for  the  various 
stages  of  testing.  First,  the  method  for  generating  dry  ice  particles  is  discussed, 
because  it  is  vital  to  this  research  effort.  Next,  the  experiments  in  particle  size  analysis 
and  PIV  conducted  at  AF1T  are  discussed.  Finally,  the  testing  conducted  at  AFRL 
using  the  TGF  is  discussed. 

3.1  Particle  Generation 

Liquid  CO2  (LCO2)  stored  in  commercially  available  dewars  is  used  throughout  the 
research.  The  dewars  measure  180  liters  in  volume  and  are  pressurized  at  approximately 
2413  kPa  at  room  temperature.  The  LC02  evaporates  into  gaseous  C02  inside  the  dewar 
which  in  turn  supplies  the  necessary  head  pressure  to  maintain  the  remaining  CO2  in 
the  liquid  form.  The  pressurized  LCO2  flows  through  a  siphon  under  its  own  pressure, 
through  a  length  of  high  pressure  flexible  stainless  steel  Swagelok®  hose  as  well  as 
stainless  steel  or  brass  compression  fittings,  attaching  to  the  feed  tube.  The  feed  tubes 
all  measure  1/8  inch  outer  diameter  (OD)  and  varied  in  ID.  The  feed  tube  was  then 
connected  to  varying  sizes  and  types  of  shroud  tubes.  Some  of  these  tubes  are  shown 
in  Figure  3.1.  This  junction  between  the  feed  and  shroud  tube  is  the  point  at  which 
the  expansion  process  ideally  begins,  allowing  the  LCO2  to  rapidly  decrease  pressure 
and  temperature,  forming  solid  CO2  particles.  While  resident  in  the  shroud  tube,  the 
particles  agglomerate  and  increase  in  size  until  they  exit  the  tube. 

Two  different  types  of  shroud  tube  were  used:  a  simple  stainless  steel  cylindrical 
tube  and  a  heat  exchanging  concentrical  tube  inside  a  tube  with  a  Kenics®  static  mixer 
element  inside  the  inner  tube.  The  static  mixer,  shown  in  Figure  3.2,  is  commonly  used 
in  the  process  industries  to  reduce  particle  size  in  liquid-liquid  and  liquid-gas  dispersions 
[19].  There  is  much  less  information  on  its  use  in  a  solid-gas  dispersion,  however.  The 
alternating  helical  mixing  elements  direct  the  flow  to  the  outer  wall  and  thoroughly 
mix  the  fluid,  as  well  as  increase  heat  transfer  to  the  outer  wall.  Temperature  controlled 
walls  through  the  use  of  a  heat  exchanging  fluid  bath  can  increase  or  decrease  this  energy 
transfer  rate,  thereby  affecting  particle  size. 
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Figure  3.1:  Feed  and  shroud  tubes  used  during  experimentation. 
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3.2  Particle  Size  Analysis 

Particle  size  data  was  collected  using  the  Malvern  Spraytec  particle  size  analyzer 
identical  to  that  used  by  Greene  and  Delapp  [7,  9].  The  Spraytec  was  set  up  in  an  open 
air  environment,  exposed  to  ambient  temperatures  and  pressures.  Different  injection 
tube  combination  were  used  to  obtain  comparative  size  measurements. 

Table  3.1  shows  the  different  combinations  of  feed  tubes  and  shroud  tubes  used, 
as  well  as  fluid  temperature  of  the  mixing  tube  heat  exchanger.  In  each  case  the  exit 
plane  of  the  shroud  tube  was  held  at  a  constant  relative  position  to  the  laser  beam.  The 
perpendicular  distance  from  the  tube  exit  to  the  laser  was  kept  at  four  inches  and  the 
intersection  of  the  injector  longitudinal  axis  and  the  laser  beam  axis  was  kept  at  5.5 
inches  from  the  optics  housing,  as  shown  in  Figure  3.3.  Greene  found  that  placing  the 
particle  stream  closer  to  the  receiver  optics  provided  more  reliable  results  [10]. 

The  particles  being  analyzed  were  solid  CO2  at  temperatures  well  below  the  am¬ 
bient  air  temperatures,  so  even  at  low  relative  humidity,  water  vapor  in  the  air  could 
condense  on  the  particles,  increasing  their  size.  To  minimize  the  effect  of  humidity  on 
the  particle  size  analysis,  a  coflow  of  conditioned  dry  air  was  supplied  through  an  Alicat 
Scientific  mass  flow  controller.  The  supply  pressure  of  the  dry  air  limited  the  coflow  to 
47  SLPM  through  the  mass  flow  controller.  This  air  supply  also  provided  clean  air  for 
the  Spraytec  optics,  reducing  the  air  available  to  minimize  humidity  effects  on  dry  ice 
crystals  through  the  laser  beam. 


Table  3.1:  Feed  and  shroud  tube  sizes 


Shroud  Tube  ID  (in) 

Feed  Tube  ID  (in) 

Temperature  (°C) 

0.03 

N/A 

0.33 

0.055 

N/A 

0.068 

N/A 

0.055 

N/A 

0.068 

N/A 

0.31  (Mixing) 

0.068 

30 

0.068 

40 

0.068 

50 

0.40 

0.068 

N/A 
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Figure  3.3:  Malvern  Spraytec  particle  size  analyzer  setup. 

Operation  of  the  Spraytec  was  relatively  straightforward.  After  system  initializa¬ 
tion,  a  background  noise  measurement  was  taken  with  the  laser  and  air  supply  on  and 
without  particles  present.  Background  noise  on  each  receiving  ring  is  subtracted  auto¬ 
matically  from  the  signal  during  actual  measurements  to  improve  accuracy.  When  ready 
to  perform  measurements,  a  time  history  hie  was  opened  and  the  laser  turned  on.  A 
C02  dewar  was  connected  by  hose  to  the  injector  tube  via  stainless  steel  and/or  brass 
compression  fittings  and  the  dewar  valve  was  opened  to  allow  how  of  liquid  C02  to  the 
tubes.  Each  test  point  was  run  between  30  seconds  and  one  minute,  depending  on  the 
observed  steadiness  of  the  how.  In  some  instances,  more  likely  with  the  mixing  tubes, 
the  solid  and  gaseous  C02  would  sporadically  clog  for  seconds  at  a  time  and  then  release, 
requiring  more  time  to  collect  sufficient  data. 

For  test  points  using  the  static  mixing  shroud  tube,  a  Thermo  Scientihc  NESLAB 
RTE  Series  refrigerated  bath  as  seen  in  Figure  3.4  supplied  fluid  how  through  the  heat 
exchanging  element  of  the  shroud  to  maintain  a  constant  shroud  temperature.  Deionized, 
distilled  water  was  fed  through  hexible  plastic  tubing  to  the  heat  exchanger  and  then 
back  to  the  bath  to  be  reconditioned  at  the  desired  temperature. 
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Figure  3.4:  Refrigerated  bath  filled  with  deionized,  distilled  water  for  heat  exchange. 

3.3  AFIT  Wind  Tunnel  Tests 

3.3.1  Wind  Tunnel.  Prior  to  its  use,  the  AFIT  variable  throat  6  inch  by  6.5 
inch  supersonic  blowdown  tunnel  shown  in  Figure  3.5  was  reconfigured  with  a  new  test 
section  to  provide  cleaner  flow  than  previously  present.  The  variable  throat  was  held 
constant  throughout  testing  and  set  so  that  Mach  1.92  flow  was  maintained  through  the 
test  section.  For  similarity  to  previous  work,  a  10°  ramp  was  chosen  to  provide  a  two- 
dimensional  oblique  shock  for  study.  Rather  than  redesign  the  tunnel  floor  to  incorporate 
this  ramp,  a  rapid  prototyping  Machine  was  used  to  build  a  ramp  that  would  fit  in  the 
existing  tunnel  floor.  Once  configured  with  a  10°  ramp,  the  tunnel  was  ready  for  use. 

The  tunnel  was  supplied  with  dry  air  by  a  tank  pressurized  to  approximately  180 
psi.  A  gate  valve  maintained  a  seal  upstream  of  the  wind  tunnel  when  not  in  use.  To 
operate  the  tunnel,  this  gate  valve  was  opened  manually.  This  pressurized  the  supply 
line  to  the  tunnel  which  was  then  kept  pressurized  by  a  pneumatic  butterfly  valve.  The 
butterfly  valve  was  operated  and  electronically  controlled  by  a  Fisher®  system.  The 
user  could  specify  a  set  point  pressure  for  the  stagnation  chamber  and  through  use  of 
a  pressure  transducer  in  the  chamber,  the  controller  would  operate  the  butterfly  valve 
to  maintain  the  set  point  pressure.  When  the  supply  pressure  dropped  below  the  spec¬ 
ified  set  point  pressure,  the  valve  would  close  automatically.  If  desired,  the  user  could 
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Figure  3.5:  AFIT  6  inch  by  6.5  inch  supersonic  blowdown  wind  tunnel. 


manually  open  and  close  the  valve  as  well.  When  collecting  PIV  data  during  this  exper¬ 
iment,  a  set  point  pressure  of  40  psi  was  used.  Once  initiated,  the  process  of  opening 
the  butterfly  valve  took  approximately  10-15  seconds  to  reach  the  set  point  pressure  and 
then  maintained  this  condition  for  approximately  20  seconds  before  the  air  supply  was 
exhausted  and  the  valve  closed. 

3.3.2  Schlieren.  Schlieren  video  was  used  to  monitor  the  upstream  flow  condi¬ 
tions  and  to  assess  whether  or  not  the  new  test  section  of  the  tunnel  provided  clean  flow. 
After  several  runs  of  the  tunnel,  however,  it  became  apparent  that  the  schlieren  provided 
the  best  means  of  determining  when  the  supersonic  flow  was  steady.  With  limited  run 
time  of  the  tunnel,  it  was  necessary  to  obstruct  the  field  of  view  (FOV)  of  the  schlieren  so 
that  it  did  not  impact  PIV  data  collection.  Therefore,  the  schlieren  setup  only  captured  a 
small  portion  of  the  flow  immediately  upstream  of  the  PIV  FOV.  The  approximate  FOV 
is  shown  in  Figure  3.6.  Light  for  the  schlieren  setup  was  emitted  from  a  mercury  short 
arc  HBO®  light  that  passed  through  a  lens  to  focus  the  light  down  to  a  point  through 
an  iris  to  pass  only  this  point  source  of  light.  The  iris  was  placed  at  the  focal  point  of 
a  concave  mirror  to  reflect  and  collimate  the  light  for  passage  across  the  wind  tunnel 
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Figure  3.6:  Approximate  location  of  the  field  of  view  for  the  schlieren  and  PIV  camera. 
Both  a  weak  shock  caused  by  a  junction  in  the  tunnel  floor  and  the  shock  caused  by  the 
wind  tunnel  ramp  are  visible. 

test  section.  The  mirrors  used  were  both  10  inches  in  diameter  with  a  focal  length  of  80 
inches.  The  second  mirror  reflected  the  light  down  to  a  point  at  which  a  vertical  knife 
edge  was  placed  to  reduce  the  light  intensity  evenly  across  the  image,  allowing  shadow 
visualization  of  the  density  gradients  in  the  light  path.  A  Photron  Fastcam-X  1280  PCI 
camera  captured  the  images  at  a  frame  rate  of  1000  frames  per  second  and  an  exposure 
duration  of  approximately  7.8  /xs,  and  then  they  were  recorded  by  the  Photron  software 
system  on  a  Compaq  computer.  Recording  was  manually  initiated  at  the  start  of  each 
run  and  halted  when  the  system  reached  its  capacity.  The  setup  of  the  schlieren  followed 
the  z-type  Herschellian  system  shown  in  Figure  2.5. 

3.3.3  PIV  System.  A  Dantec  Dynamics  PIV  system  was  used  to  acquire  data 
in  the  AFIT  wind  tunnel.  The  laser  source  was  a  a  New  Wave  Research  Solo  200XT 
Nd: YAG  frequency  doubled  laser  operating  at  532  nm  and  a  maximum  power  of  200mJ. 
An  optical  arm  consisting  of  a  series  of  optics  was  used  to  manipulate  the  laser  beam 
for  positioning  above  the  wind  tunnel  and  to  spread  the  beam  into  a  planar  laser  sheet 
approximately  4mm  wide.  Images  were  captured  by  a  FlowSense  4M  CCD  camera  with 
a  2048  x  2048  pixel  array  outfitted  with  a  Nikkor  60  mm  f2.8  lens  and  then  fed  to 
the  Dantec  Dynamic  Studio  software  package  for  recording.  The  software  also  provided 
control  of  the  camera  and  laser  triggering  functions.  This  system  was  set  up  to  capture 
particle  motion  in  air  flow  across  a  shock  generated  by  a  10°  ramp  using  a  vertical  laser 
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sheet  aligned  with  the  focal  plane  of  the  camera,  positioned  perpendicular  to  the  flow 
direction,  with  flow  from  left  to  right  as  seen  in  Fig.  3.7.  After  positioning  the  laser  and 
optical  arm  such  that  the  laser  sheet  was  located  in  the  center  of  the  tunnel  and  aligned 
with  the  flow  direction,  a  calibration  grid  was  aligned  with  the  laser  sheet  to  focus  the 
camera  and  scale  the  held  of  view.  Figure  3.8  shows  the  calibration  image  used.  The 
grid  imaged  is  of  known  distance  so  a  scale  factor  can  be  calculated  to  determine  the 
relative  size  each  pixel  captures  on  the  focal  plane.  Performing  this  process  resulted  in 
an  image  of  the  flow  held  with  a  held  of  view  that  measured  77.8  mm  by  77.8  mm,  giving 
a  resolution  of  26.3  pixels  per  millimeter.  This  translated  to  a  scale  factor  of  5.096  in 
the  Dantec  software.  To  ensure  the  camera  was  properly  focused  and  to  test  the  system 
setup,  an  injector  was  placed  on  top  of  the  wind  tunnel  and  angled  downward  to  direct 
the  CO2  particles  into  the  test  section  through  the  PIV  FOV.  Figure  3.9  shows  this 
arrangement.  For  wind  tunnel  runs  and  PIV  data  collection,  image  pairs  were  collected 
using  a  1  /is  delay  and  a  repetition  frequency  between  2-4  Hz.  For  the  how  conditions  of 
the  wind  tunnel,  the  system  timing  equated  to  a  particle  displacement  of  about  12  pixels 
in  the  freestream. 
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Figure  3.7:  Laser  and  camera  setup  for  PIV  measurements  in  the  AFIT  6  inch  by  6.5 
inch  wind  tunnel. 
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Figure  3.8:  Calibration  grid  imaged  to  scale  the  field  of  view. 
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Figure  3.9: 


Injection  from  above  the  test  section  used  for  system  checkout. 
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The  recipe  for  processing  the  image  pairs  began  with  a  cross-correlation  using  two 
different  interrogation  region  (IR)  sizes:  64  by  64  pixels  (2.43  mm  by  2.43  mm)  and  32  by 
32  pixels  (1.22  mm  by  1.22  mm).  The  full  image  was  processed  with  these  IR  sizes  using 
a  50%  overlap  for  the  322  pixel  IR  and  a  75%  overlap  for  the  642  pixel  IR.  A  Gaussian 
window  function  using  equation  3.1  was  applied  to  attenuate  the  signal  strengths  closer  to 
the  edge  of  the  IR,  thereby  weighting  the  correlation  towards  the  signal  near  the  center 
of  the  IR  [14].  In  equation  3.1  the  coefficient  k  determines  the  width  of  the  window 
function,  and  for  this  research  k  =  1.19  was  used.  The  weight,  W,  was  determined  using 
the  pixel  location,  (■ m,n ),  and  the  IR  size,  ( M,N ).  A  No-DC  filter  was  also  applied  to 
remove  the  background  noise  to  effectively  enhance  the  signal-to-noise  ratio  [14].  The 
correlation  was  followed  by  a  peak  validation  using  a  value  of  1.1  and  a  range  validation 
to  remove  outlier  vectors.  Since  the  range  of  flow  speeds  were  known  from  pressure 
transducer  data  in  the  tunnel  as  well  as  visual  correlation  of  particles  in  the  image  pairs, 
the  range  of  vectors  passed  was  restricted  to  200-600  m/s  for  the  streamwise  component 
of  velocity  (U)  and  -100  to  100  m/s  for  the  vertical  component  of  velocity  (V).  Finally, 
averages  were  calculated  for  each  condition  using  only  vectors  that  had  passed  both  the 
peak  and  range  validation  methods. 


W  ( m ,  n)  =  exp 


(3.1) 


In  equation  3.1,  W  =  weight  value;  (m,n)  =  pixel  location;  k  =  window  function  width 
coefficient;  and  (M,N)  =  interrogation  region  size. 


3-4  AFRL  Wind  Tunnel  Tests 

3. 4-1  Trisonic  Gas-dynamics  Facility.  The  AFRL  Trisonic  Gas-dynamics  Fa¬ 
cility  (TGF)  is  a  closed  loop,  variable  density,  continuous  flow  wind  tunnel  with  a  2  feet 
by  2  feet  (0.61  m  by  0.61  m)  test  section  (Figure  3.10).  It  uses  interchangeable  nozzle 
blocks  to  operate  in  the  subsonic,  transonic  and  supersonic  regimes.  In  this  test  only 
the  subsonic  nozzle  block  was  used,  providing  the  performance  characteristics  seen  in 
Figure  3.11.  The  tunnel  was  operated  by  Shawn  Raisch,  Paul  Olekas,  and  Rick  Zimmer- 


man  from  Jacobs  Engineering.  A  LaVision  PIV  system  was  arranged  and  operated  by 
Dr.  Ryan  Schrnit  with  assistance  from  Tom  Presdorf,  PE  (AFRL/RB)  and  Dr.  Aaron 
Altman  (University  of  Dayton).  Dry  ice  experiments  were  performed  on  the  same  days 
as  tests  of  a  Vicount  system,  beginning  on  29  January  2010  and  ending  on  5  February 
2010.  The  tunnel  was  conhgured  with  CO2  injectors  in  four  locations;  three  in  the 
stagnation  chamber  and  one  just  upstream  of  the  test  section.  Figure  3.12  shows  two 
of  the  three  injecting  into  the  wall  of  the  stagnation  chamber  and  Figure  3.13  is  a  view 
looking  upstream  from  the  test  section  showing  the  fourth  injector.  The  TGF  was  used 
on  two  separate  days,  first  to  assess  successful  generation  and  persistence  of  particles 
through  the  test  section,  and  second  to  perform  PIV  using  the  CO2.  For  day  one,  three 
combinations  of  feed  and  shroud  tubes  were  used,  each  at  a  different  location  in  the 
tunnel.  Each  injector  was  used  individually  at  each  test  condition  shown  in  Table  3.2. 
The  best  seeding  observed  during  this  test  would  then  be  used  for  PIV  on  the  second 
day. 


Table  3.2:  Test  conditions  for  TGF. 


Mach  Number 

Pressure  (PSF) 

0.89 

255 

0.79 

501 

500 

0.5  (approx.) 

1000 

1500 

0.30 

516 

0.30 

1006 

0.30 

1506 

3. 4-2  PIV  System.  The  PIV  system  used  for  TGF  testing  consisted  of  a 
LaVision  computer  and  software  package,  a  PCO  camera,  and  a  New  Wave  Research 
Solo  PIV  Nd:YAG  laser.  For  the  first  day  the  system  was  set  up  with  two  cameras,  one 
on  either  side  of  the  tunnel,  capturing  different  regions  of  the  tunnel.  Because  it  was  not 
known  where  the  particles  would  appear  in  the  tunnel  held  of  view,  the  laser  provided  a 
sheet  of  light  across  the  tunnel  vertically  such  that  the  how  direction  was  normal  to  the 
laser  sheet.  This  provided  the  maximum  coverage  of  the  tunnel  test  section  to  strictly 
observe  particles,  thus  no  velocimetry  data  was  collected.  For  the  second  day  of  testing, 
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Figure  3.11:  TGF  subsonic  test  section  performance  chart. 

[20] 

the  laser  sheet  was  rotated  to  be  in  plane  with  the  flow  direction.  The  geometry  of  the 
TGF  limited  the  camera  FOV  such  that  the  only  way  to  collect  images  perpendicular 
to  the  laser  sheet  was  to  bring  the  laser  sheet  through  the  tunnel  at  an  angle.  In 
this  arrangement,  the  camera  was  approximately  90°  from  the  laser  plane,  though  this 
plane  did  not  align  with  standard  coordinate  frames.  At  each  test  condition  shown  in 
Table  3.2,  the  LaVision  system  captured  100  image  pairs.  The  image  pairs  were  then 
processed  using  a  similar  cross-correlation  technique.  Differences  in  the  software  did  not 
allow  identical  post  processing,  but  in  this  case  a  range  validation  was  applied  to  the 
streamwise  velocity  component.  Vectors  with  a  streamwise  velocity  outside  the  ranges 
specified  in  Table  3.3  below  were  not  used  in  the  mean  velocity  computations. 

Table  3.3:  Range  validation  for  TGF  data. 

Min  (m/s)  Expected  Mach/speed  Max  (m/s) 

25  0.3/100  175 

100  0.5/170  300 

150  0.8/270  350 


Figure  3.12:  Exterior  view  of  CO2  injectors  in  the  stagnation  chamber  wall  of  the  TGF 
with  an  interior  view  inset  looking  up  from  the  access  door. 
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IV.  Results  and  Discussion 


Results  for  testing  conducted  at  AFIT  and  AFRL  are  presented  below.  First, 
the  particle  size  analysis  conducted  at  AFIT  is  presented.  Based  on  these  re¬ 
sults,  particular  configurations  were  chosen  for  insertion  into  the  AFIT  wind  tunnel. 
PIV  results  are  next  presented  for  the  AFIT  tunnel  testing.  Finally,  results  for  testing 
conducted  at  the  AFRL/RB  TGF  is  presented. 

4-1  Particle  Size  Analysis 

Particle  size  analysis  was  accomplished  using  the  Malvern  Spraytec.  The  CO2  was 
run  approximately  45  seconds  through  the  laser  for  different  configurations  of  feed  and 
shroud  tubes  and  data  was  recorded  by  the  RTSizer  software.  Because  of  the  differences 
in  environments  between  where  the  Spraytec  was  used  and  where  the  particles  will  be 
needed  (wind  tunnel),  the  particle  size  information  presented  is  for  relative  compara¬ 
tive  information  rather  than  for  true  particle  size  information.  In  each  case  the  Sauter 
mean  diameter,  D32,  was  chosen  to  give  representative  particle  size  information  based 
on  Greene’s  results  that  showed  this  was  a  good  indicator  of  size  [10].  To  arrive  at 
each  value,  random  snapshots  of  the  recorded  data  were  exported  and  the  mean  and 
standard  deviation  of  these  discrete  measurements  were  calculated.  Particle  dispersion 
is  addressed  and  the  effects  due  to  the  feed  tube  ID,  shroud  tube  ID,  shroud  tube  type, 
and  heat  exchanging  fluid  bath  temperature  are  compared  below. 

Though  mean  diameter  data  is  presented,  it  is  not  without  error  in  addition  to 
the  standard  deviation  of  the  data  points  shown.  The  CO2  particle  size  distribution  was 
nonuniform,  so  even  though  a  single  number  is  presented  to  describe  the  particles  at  each 
test  condition,  the  actual  particles  varied  significantly.  In  most  cases,  the  percent  number 
of  particles  was  greater  for  the  smaller  particles  than  larger  ones,  but  larger  particles 
made  up  the  majority  of  the  percent  volume  measurement.  The  percent  number  and 
percent  volume  measurements  from  the  size  analyzer  were  used  to  determine  the  Sauter 
mean  diameter,  but  it  is  just  a  calculation  of  a  particularly  sized  particle  spray  that 
would  have  the  same  volume  to  surface  area  ratio  as  the  actual  spray.  Figures  4.1  and 
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4.2  are  typical  distributions  of  the  particle  number  and  volume.  So,  even  though  the 
mean  particle  size  was  in  the  range  presented,  a  large  variation  of  size  was  present. 


Figure  4.1:  Percent  of  total  number  of  particles  measured  at  each  particle  diameter. 


4-1.1  Particle  Dispersion.  Two  different  types  of  shroud  tube  were  used  for 
measurements:  a  simple  shroud  tube  (Figure  4.3)  and  a  static  mixing  shroud  tube  (Figure 
4.4).  There  was  a  strong  difference  in  the  dispersion  of  CO2  as  it  exited  the  tube.  Figures 
4.3  and  4.4  show  the  visible  difference  in  the  spread  angle  as  the  particles  sprayed  out. 
In  Figure  4.3  the  particles  and  gaseous  CO2  were  sprayed  predominantly  in  the  axial 
direction  whereas  in  Figure  4.4  the  particles  and  gaseous  CO2  exited  at  a  slower  velocity 
and  spread  much  more  rapidly. 

4-1.2  Effect  due  to  Feed  Tube  Size.  Three  different  diameters  of  feed  tubes 
were  used:  0.03  inch,  0.055  inch,  and  0.068  inch.  All  feed  tubes  for  size  analysis  were 
12.1  inches  in  length.  Each  tube  was  married  with  the  0.33  inch  ID  simple  shroud 
tube  to  isolate  the  effects  of  the  feed  tubes.  Figure  4.5  shows  the  change  in  the  Sauter 
mean  diameter,  D32,  for  the  three  tubes  with  error  bars  at  one  standard  deviation.  It 
is  apparent  that  an  increase  in  the  ID  size  of  the  feed  tube  decreases  the  mean  particle 
size  significantly,  which  is  in  agreement  with  the  results  of  Greene.  An  effect  of  changing 
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Figure  4.2: 


Percent  of  total  volume  particles  of  a  given  size  make  up. 
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Figure  4.4:  Characteristic  profile  of  particles  exiting  a  mixing  shroud  tube. 

orifice  size  that  was  not  measured  but  must  be  mentioned  is  the  increase  in  mass  flow  rate 
of  the  CO2  with  increasing  orifice  size.  A  larger  amount  of  liquid  CO2  flowed  through  the 
orifice  and  expanded,  forming  more  particles,  but  also  caused  the  velocity  in  the  shroud 
to  increase.  Conversely,  a  small  ID  feed  tube  restricted  the  mass  flow,  thus  slowing  the 
flow  through  the  shroud  tube.  As  described  above,  a  decrease  in  the  speed  through  the 
shroud  tube  allowed  for  larger  particles  to  grow,  explaining  the  larger  sizes  seen  with  the 
0.030  inch  ID  feed  tube. 

4-1-3  Effect  due  to  Shroud  Size.  Two  different  diameters  of  shroud  tube  ID 
were  used:  0.33  inch  and  0.40  inch.  Each  of  these  tubes  measured  10  inches  in  length. 
Figure  4.6  gives  the  mean  D32  and  standard  deviation  for  these  two  configurations  using 
the  0.068  inch  feed  tube.  With  a  slightly  larger  diameter  and  a  constant  mass  flow  rate 
as  governed  by  the  orifice,  the  residence  time  of  particles  was  higher  for  the  larger  shroud 
tube.  This  result  is  also  consistent  with  those  reported  by  Greene. 

4-1-4  Effect  due  to  Shroud  Type.  In  addition  to  considering  the  effect  of  shroud 
diameter,  two  different  types  of  shroud  tube  were  studied:  the  simple  straight  tubes 
used  above  and  the  static  mixing  tube.  The  mixing  tube  had  slightly  thicker  walls  than 
the  similar  sized  simple  tube,  which  reduced  the  ID  from  0.33  inch  for  the  simple  tube 
to  0.31  inch  for  the  mixing  tube.  Also,  the  tube  length  was  different;  the  simple  tubes 
were  10  inches  long  whereas  the  mixing  tubes  were  14  inches  long.  Greene  found  that 
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Figure  4.5: 


Figure  4.6: 


Average  D[3 ,2]  for  Varying  Feed  Tube  Inner  Diameters 


Feed  tube  Inner  Diameter  (inches) 


D32  for  varying  feed  tube  fD  and  a  constant  0.33  inch  fD  shroud  tube. 


Average  D[3][2]  Using  0.068"  ID  Feed  Tube  and  Varying  Shroud  Size 
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D32  for  varying  shroud  tube  size  and  constant  0.068  inch  feed  tube  ID. 
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particle  size  plateaued  at  a  length  of  approximately  10  inches  and  did  not  vary  largely 
for  a  shroud  tube  20  inches  long  so  this  difference  was  assumed  to  have  little  effect  [10]. 
The  main  differences  were  two  fold:  the  addition  of  the  Kenics®  static  mixer  as  well 
as  the  heat  exchanging  fluid  bath.  With  no  fluid  present,  Figure  4.7  shows  the  mean 
D32  and  standard  deviation  for  the  two  types  of  tubes.  These  results  suggest  that  the 
Kenics®  mixer  performs  as  expected  in  reducing  the  mean  particle  size,  D32,  by  a  factor 
of  approximately  three  even  with  no  fluid  bath  present. 
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Average  D[3][2]  Using  0.055"  ID  Feed  Tube  and  Varying  Shroud  Type 


Mixing  Shroud 


Shroud  Type 


Straight  Shroud 


Figure  4.7:  D32  for  a  simple  versus  mixing  shroud  tube  type. 


4-1-5  Effect  due  to  Thermal  Bath  Temperature.  Using  deionized  and  distilled 
water,  a  refrigerated  bath  heated  or  cooled  the  water  to  a  designated  temperature  and 
pumped  the  water  through  the  heat  exchanging  shroud  during  data  acquisition.  Each 
data  acquisition  event  lasted  approximately  one  minute.  In  that  time,  the  fluid  bath 
temperature  remained  nearly  constant,  dropping  less  than  0.5°C.  However,  at  the  lower 
temperature  measurements  (20°C  and  30°C),  the  water  froze  in  the  shroud  housing  such 
that  the  temperature  of  the  heat  exchanging  fluid  was  not  maintained  in  the  shroud  at 
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the  desired  setting.  Figure  4.8  shows  the  results  of  the  different  bath  temperatures  using 
the  0.068  inch  feed  tube.  It  was  evident  that  the  small  temperature  range  of  20°C  -  50°C, 
as  operated,  produced  no  significant  trends  in  particle  size. 


Average  D[3.2]  Using  0.068"  ID  Feed  Tube  and  0  31"  ID  Mixing  Shroud 
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Figure  4.8:  D32  for  the  mixing  shroud,  0.068  inch  ID  feed  tube  and  varying  shroud 

bath  temperature. 
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30  40 
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4-1-6  Sources  of  Error  in  Particle  Size  Analysis.  As  discussed  at  the  beginning 
of  this  section,  the  Sauter  mean  diameter  data  presented  is  not  the  actual  particle  size, 
but  rather  an  equivalent  sized  sphere  that  represents  the  particle  size  distribution.  It  is 
affected  by  the  limitations  of  the  particle  size  analyzer,  being  limited  to  the  detection  and 
calculation  of  a  range  of  particles.  The  manufacturer  suggested  measurements  less  than 
about  2  /irn  are  not  reliable.  As  is  evident  in  Figure  4.1,  the  Spraytec  system  measured 
the  majority  of  particles  in  a  range  less  than  this,  so  there  is  uncertainty  in  this  measure. 
Also,  the  experimental  setup  is  an  error  source  in  that  it  is  conducted  in  an  environment 
that  is  exposed  to  ambient  air,  allowing  humidity  to  affect  the  particle  size.  Water  vapor 
can  condense  on  particles,  increasing  their  size  during  measurements.  Conditioned  dry 
air  was  provided  in  a  coflow  to  reduce  this  effect,  however,  it  may  have  still  been  present. 
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Finally,  the  flow  of  C02  was  not  steady.  Clogging  reduced  the  flow  significantly  at  times, 
causing  a  decrease  in  particles.  This  issue  was  more  prominent  with  the  mixing  tube, 
requiring  more  time  to  collect  size  measurements. 

4-2  AFIT  Wind  Tunnel  Results 

The  AFIT  6  inch  by  6.5  inch  supersonic  blowdown  tunnel  was  configured  with  a 
10°  ramp  to  observe  the  effects  across  an  oblique  shock.  Temperature  and  pressure  trans¬ 
ducers  installed  on  the  tunnel  were  used  to  collect  corresponding  data  for  the  tunnel  flow 
conditions.  This  data  was  then  used  to  calculate  the  speed  of  the  flow  and  approximate 
the  shock  wave  angle  using  equation  2.11  and  isentropic  flow  theory.  Using  a  set  point 
of  40  psi,  the  approximate  tunnel  speed  was  calculated  to  be  Mach  1.92  or  475  m/s.  At 
this  speed  the  10°  ramp  would  produce  a  41°  shock  wave.  These  values  were  used  for 
comparison  to  PIV  data  collected.  Both  temperature  and  pressure  changed  during  the 
course  of  each  tunnel  operation  so  there  was  some  variability  in  the  speed  and  shock 
angle  that  was  unavoidable.  Figure  4.9  shows  the  typical  change  in  conditions  for  each 
operation  of  the  tunnel. 

With  the  time  delay  between  laser  pulses  set  to  1  /xs  and  with  64  pixels  correspond¬ 
ing  to  2.43  mm,  particles  in  the  freestreeam  typically  traveled  about  12  pixels.  Most  PIV 
systems  are  operated  such  that  particles  travel  no  more  than  about  0.25  times  the  length 
of  the  IR  [3,  21].  Thus,  the  64  by  64  pixel  subregion  size  was  expected  to  be  reasonable, 
and  a  32  by  32  subregion  size  was  also  explored. 

4-2.1  Simple  Shroud  Tubes.  The  tunnel  was  configured  with  three  dewars  and 
three  injection  points,  each  using  a  0.055  inch  ID  feed  tube  and  0.33  inch  simple  shroud 
tube.  Laser  power  was  set  at  100  mJ  per  pulse  and  the  camera  aperture  was  set  at  f5.6. 
Unfortunately  due  to  tunnel  vibrations  and  the  duration  of  operation  at  relatively  stable 
conditions,  only  10-20  image  pairs  were  collected  during  each  run.  Several  tunnel  runs 
were  combined  to  accumulate  174  image  pairs.  Using  the  process  described  in  chapter 
three,  the  image  pairs  were  correlated  to  produced  vector  statistics  across  the  FOV. 
Figures  4.10  and  4.11  show  the  horizontal  (U)  velocity  component  using  an  IR  size  of  32 
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Figure  4.9:  Typical  transient  conditions  of  AFIT  blowdown  wind  tunnel  operation. 
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Temperature  (K) 


by  32  pixels  and  64  by  64  pixels,  respectively.  While  the  results  are  generally  similar, 
the  64  by  64  IR  size  led  to  a  smoother  plot.  The  freestream  velocity  in  both  cases 
corresponded  well  with  the  calculated  value  of  475  m/s.  The  velocity  downstream  of  the 
shock  also  correlated  well  with  a  calculated  value  of  439  m/s. 


20  40  60  80 
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Figure  4.10:  Horizontal  velocity  component  using  simple  tubes  and  a  32  x  32  pixel  IR. 

The  vertical  velocity  (V)  perpendicular  to  the  flow  direction  is  shown  for  the  two 
different  IR  sizes  in  Figures  4.12  and  4.13.  In  this  case,  the  V  component  should  be  zero 
upstream  of  the  shock  but  the  PIV  data  apparent  in  both  figures  showed  a  slight  positive 
vertical  component  of  between  10-20  m/s.  This  equated  to  approximately  1-2  degrees 
of  angularity  in  the  flow  upstream  of  the  shock.  Downstream  of  the  oblique  shock  the 
calculated  velocity  from  shock  relations  is  77  m/s  while  the  PIV  data  suggested  a  value 
closer  to  60  m/s.  This  may  be  related  to  the  apparent  flow  angularity.  By  observing  the 
schlicren  images  it  was  apparent  that  the  flow  upstream  of  the  shockwave  caused  by  the 
ramp  was  not  entirely  free  of  disturbances.  It  was  possible  that  these  disturbances  led 
to  the  observed  angularity  of  the  flow.  Figure  4.14  shows  a  weak  oblique  shock  or  Mach 
wave  upstream  of  the  shock  caused  by  the  ramp.  Another  factor  that  could  account  for 
this  was  that  this  vertical  motion  may  have  been  a  result  of  the  minimum  resolution  of 
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Figure  4.11:  Horizontal  velocity  component  using  simple  tubes  and  a  64  x  64  pixel  IR. 


the  PIV  processing.  For  the  PIV  FOV  and  the  camera  resolution,  a  motion  of  one  pixel 
in  an  image  pair  equated  to  about  38  m/s. 
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Figure  4.12: 


Figure  4.13: 


Vertical  velocity  component  using  simple  tubes  and  a  32  x  32  pixel  IR. 
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Vertical  velocity  component  using  simple  tubes  and  a  64  x  64  pixel  IR. 
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Weak  Shock 


Ramp  Induced 
Shock 


Figure  4.14:  Schlieren  image  of  flow  upstream  of  PIV  FOV  showing  a  weak  shock 

caused  by  tunnel  roughness  and  a  strong  shock  caused  by  the  10°  ramp. 
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Having  the  horizontal  and  vertical  velocity  components  allows  the  calculation  of 
the  velocity  normal  to  the  shock,  Vnorm,  using  equation  4.1  for  each  IR. 

Vfiorm  —  U  sin  V  cos  /3  (4.1) 

Figures  4.15  and  4.16  show  the  results  of  the  32  and  64  pixel  IR  correlations  respectively. 
The  velocity  tracked  well  with  the  calculated  upstream  and  downstream  velocities  normal 
to  the  shock:  312  m/s  and  230  m/s.  By  selecting  vectors  in  a  direction  normal  to 
the  shock,  the  particle  response  could  be  compared  with  the  theoretical  response  as 
determined  using  Melling’s  method  described  in  chapter  2.  In  Figure  4.17  the  response 
of  the  particles  measured  using  both  IR  sizes  is  compared  with  the  theoretical  response 
of  different  size  particles.  The  PIV  data  for  both  followed  the  trend  of  a  2.0  jum  particle 
with  a  density  of  1.18  g/cm3,  and  the  smaller  IR  size  did  not  show  any  advantage  when 
determining  this.  Using  the  smaller  IR  size  could  have  improved  the  resolution  and 
altered  the  response  curve,  but  in  this  case  it  did  not. 
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Figure  4.15:  Velocity  normal  to  the  shock  using  simple  tubes  and  a  32  x  32  pixel  IR. 


47 


20  40  60  80 

x(mm) 


Figure  4.16:  Velocity  normal  to  the  shock  using  simple  tubes  and  a  64  x  64  pixel  IR. 


Particle  Response  to  an  Oblique  Shock  using  Straight  Shroud  Tubes 


Figure  4.17:  Particle  response  using  two  IR  sizes  compared  to  theoretical  response  for 
varying  particle  sizes. 
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4-2.2  Static  Mixing  Shroud  Tubes.  For  the  second  set  of  tests  the  simple 
shroud  tubes  were  replaced  by  three  static  mixing  shroud  tubes  holding  everything  else 
constant.  Initial  runs  of  the  tunnel  in  this  configuration  resulted  with  limited  particle 
coverage  in  the  image  pairs  collected.  The  hypothesis  was  that  the  particles  were  too 
small  to  be  imaged  successfully  using  the  same  camera  and  laser  settings.  Therefore, 
the  camera  f-stop  was  increased  one  stop  to  f4.0  and  the  laser  power  was  increased  to 
200  mJ  per  pulse.  These  changes  improved  the  visibility  of  particles  in  the  images  so  all 
data  collected  on  the  mixing  shroud  tubes  used  these  settings  for  the  camera  and  laser. 
Overall,  150  image  pairs  were  collected  and  used  to  produce  the  results  below.  Figures 
4.18  and  4.19  show  the  horizontal  (U)  velocity  component  using  an  IR  size  of  32  by  32 
pixels  and  64  by  64  pixels,  respectively.  Again,  the  particles  tracked  the  flow  accurately 
far  upstream  and  far  downstream  of  the  oblique  shock  wave.  The  same  is  evident  in 
the  vertical  velocity  component  shown  in  Figures  4.20  and  4.21.  The  normal  velocity 
is  calculated  in  the  same  manner  as  above  and  is  plotted  in  Figures  4.22  and  4.23.  At 
the  higher  laser  intensity  and  with  an  increased  aperture,  the  particles  generated  by  the 
static  mixing  tubes  were  able  to  track  the  freestream  and  downstream  flow  velocities 
accurately  when  compared  to  isentropic  and  shock  relations.  Across  the  shock,  the 
particles  lag  the  fluid  velocity  similar  to  a  theoretical  2.0  /jrn  particle,  as  can  be  seen 
in  Figure  4.24.  Surprisingly,  even  though  laser  power  and  aperture  changes  suggested 
smaller  particles,  the  response  characteristics  indicated  no  significant  change  in  the  size 
of  particles  produced.  For  example,  Figure  4.25  compares  a  128  by  128  pixel  region 
collected  from  a  typical  PIV  image  using  the  simple  injectors  to  a  similar  region  using 
the  mixing  injectors.  The  particles  appear  roughly  the  same  size,  however,  the  aperture 
for  the  mixing  injectors  was  opened  to  f4.0  and  the  laser  power  was  increased  to  capture 
this  image.  This  suggested  that  there  were  either  fewer  particles  or  smaller  particles 
with  the  mixing  shroud  tubes.  Figure  4.26  shows,  however,  that  the  oblique  shock  wave 
was  not  entirely  stationary  due  to  turbulent  structures  in  the  boundary  layer  or  possibly 
due  to  acoustic  interactions  in  the  tunnel.  This  can  be  seen  with  reference  to  the  red 
horizontal  reference  line,  which  is  drawn  near  the  top  of  the  oblique  shock,  indicated  by 
the  dark  diagonal  line.  The  nine  sequential  frames  show  that  there  is  shock  motion  as 
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the  dark  diagonal  line  crosses  above  and  below  the  stationary  red  reference  line.  Shock 
foot  motion  has  been  described  by  Beresh,  et.  ah,  albeit  for  a  larger  turning  angle  for 
a  Mach  5  freestream  [22],  Similar  instability  of  a  shock  foot  has  been  reported  in  the 
divergent  section  of  nozzles  as  well  [23].  Visual  inspection  of  schlicren  images  suggest 
the  shock  motion  is  approximately  4  mm  normal  to  the  oblique  shock  wave.  Because 
PIV  results  are  averages  over  the  collection  time,  the  movement  of  the  shock  would  be 
averaged  out,  resulting  in  a  larger  particle  response  distance  based  on  PIV  results. 
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Figure  4.18:  Horizontal  velocity  component  using  mixing  tubes  and  a  32  x  32  pixel 

IR, 
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Figure  4.19: 
IR. 


Figure  4.20: 
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Horizontal  velocity  component  using  mixing  tubes  and  a  64  x  64  pixel 


Vertical  velocity  component  using  mixing  tubes  and  a  32  x  32  pixel  IR. 
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Figure  4.21: 


Figure  4.22: 
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Vertical  velocity  component  using  mixing  tubes  and  a  64  x  64  pixel  IR. 
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Velocity  normal  to  the  shock  using  mixing  tubes  and  a  32  x  32  pixel  IR. 
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Figure  4.23:  Velocity  normal  to  the  shock  using  mixing  tubes  and  a  64  x  64  pixel  IR. 


Particle  Response  to  an  Oblique  Shock  using  Mixing  Shroud  Tube 


Figure  4.24:  Particle  response  using  two  IR  sizes  compared  to  theoretical  response  for 
varying  particle  sizes 
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Figure  4.25:  128  by  128  pixel  regions  for  simple  injectors  at  f5.6  (left)  and  mixing 

injectors  at  f4.0  (right). 


Figure  4.26:  Nine  sequential  frames  using  a  red  reference  line  to  indicate  shock  dis¬ 

placement  over  time. 
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Though  150  image  pairs  were  collected  and  used  for  these  calculations,  many  of 
the  image  pairs  collected  showed  sparse  particle  coverage  when  using  the  mixing  tubes.  1 
believe  this  is  the  effect  of  CO2  accumulation  either  in  the  line  or  in  the  mixing  tube  itself 
to  the  point  where  all  flow  was  restricted.  This  clogging  was  observed  when  performing 
Malvern  particle  size  analysis  as  well.  For  Malvern  testing,  most  clogs  did  not  seem  to 
persist  more  than  10  seconds,  but  with  the  limited  duration  of  PIV  image  collection  in 
the  tunnel,  a  clog  in  a  single  injector  could  reduce  the  particle  count  by  one  third  for  all 
the  images  collected  on  that  run.  If  more  than  one  line  clogged,  obviously  the  quantity 
of  particles  would  be  further  reduced. 

In  both  Figures  4.17  and  4.24,  the  measured  velocity  using  the  smaller  IR  was 
lower  than  that  found  using  the  larger  IR.  This  suggests  using  the  smaller  IR  biases  the 
vectors  produced  towards  zero.  Use  of  the  Gaussian  filter  may  be  the  root  cause.  For 
the  1  /x s  time  delay  between  image  frames,  the  pixel  displacement  is  approximately  12 
pixels.  For  a  32  pixel  IR,  the  displacement  is  one  third  the  length  of  the  IR  and  the 
Dantec  software  recommends  less  than  this  when  using  the  Gaussian  filter  [14]. 

4-2.3  Sources  of  Error  in  AFIT  Wind  Tunnel  Testing.  Though  some  discus¬ 
sion  has  been  covered  above  on  the  potential  reasons  for  the  disagreements  between  PIV 
data  and  the  expected  velocities  based  on  tunnel  instrumentation,  several  other  areas 
can  affect  the  accuracy  of  these  comparisons.  First,  the  expected  velocities  were  based 
on  tunnel  instrumentation  that  lacked  recent,  valid  calibration.  Both  pressure  and  tem¬ 
perature  transducers  were  used  in  the  tunnel  to  provide  the  necessary  information  to 
compute  expected  conditions  upstream  of  the  shock,  but  the  only  calibration  equipment 
available  at  the  time  of  testing  required  calibration  itself.  These  transducers  had  been 
used  recently  by  Freeborn  and  Greene  for  their  experiments,  and  the  calibration  per¬ 
formed  during  this  experimentation  agreed  well  with  their  data,  but  it  is  still  a  source  of 
uncertainty  [9,  11]. 

Next,  the  modifications  to  the  tunnel  test  section  prior  to  testing  may  have  in¬ 
troduced  a  weak  shock  as  seen  in  Figure  4.14.  However,  this  apparent  shock  may  have 
been  isolated  to  one  of  the  sides  of  the  tunnel  rather  than  across  the  entire  area.  The 
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schlieren  system  used  was  not  focused,  rather,  captured  all  density  gradients  between  the 
two  mirrors  and  collapsed  them  onto  a  two  dimensional  image.  If  the  shock  was  located 
directly  upstream  of  the  focal  plane  of  the  PIV  system,  the  disagreement  between  the 
velocities  would  be  increased. 

Finally,  the  PIV  system  may  introduce  errors  that  aren’t  known.  The  velocity 
results  are  based  on  a  user-defined  calibration  correlating  a  length  scale  to  a  pixel  scale. 
The  velocity  calculations  require  sub-pixel  analysis  to  resolve  velocity  differences  less  than 
about  38  m/s  in  this  setup,  which  is  the  equivalent  speed  of  a  one-pixel  displacement  for 
the  1  /is  time  delay.  If  the  calibration  is  slightly  off,  velocity  results  could  be  biased. 

4-3  AFRL  Wind  Tunnel  Results 

The  TGF  was  used  with  a  LaVision  PIV  system  and  multiple  injection  ports  and 
types  to  demonstrate  the  usability  of  this  seeding  method.  The  TGF  uses  interchangeable 
nozzle  blocks  to  produce  supersonic  flow,  but  for  the  testing  performed,  only  the  subsonic 
nozzle  was  used,  allowing  speeds  up  to  approximately  M  =  0.8  at  variable  pressures. 
Figure  3.11  shows  the  subsonic  conditions  achievable.  The  test  conditions  used  for  testing 
are  shown  in  Table  3.2. 

4-3.1  Cross  Plane.  Figures  3.12  and  3.13  show  the  initial  setup  using  two 
injection  points  in  the  stilling  chamber  through  the  wall  and  a  third  injection  point 
just  upstream  of  the  test  section.  Table  4.1  lists  the  position  and  configuration  of  each 
injector.  To  improve  the  chances  of  seeing  particles,  each  injection  point  used  a  different 
combination  of  tubes,  introducing  two  that  had  not  been  used  for  data  collection  in  the 
AFIT  tunnel:  a  curved  0.4  inch  ID  shroud  tube  and  a  0.17  inch  ID  shroud  tube.  The 
third  injector  was  identical  to  that  used  in  the  AFIT  wind  tunnel:  0.055  inch  ID  feed 
tube  and  0.33  inch  ID  shroud  tube.  This  initial  test  was  used  to  determine  whether 
or  not  seeding  was  present  through  the  test  section,  and  if  so,  where.  The  increased 
distance  from  the  injector  port  to  the  test  section  (approximately  20  feet)  and  the  lower 
tunnel  speed  meant  that  the  particles  had  to  persist  much  longer  in  the  TGF  than  in  the 
smaller  AFIT  wind  tunnels.  Therefore,  the  PIV  system  projected  the  laser  sheet  such 
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that  the  flow  direction  was  normal  to  the  laser  plane.  Two  cameras  were  nsed  to  image 
the  left  and  right  sides  of  the  tunnel  using  Scheimphlug  mounts  to  maintain  focus  from  an 
off-axis  direction.  This  arrangement  is  shown  in  Figure  4.27.  The  images  collected  were 
not  processed  for  PIV,  rather  they  were  used  to  determine  whether  particles  persisted 
and  where  they  were  present. 


Table  4.1:  Injector  configurations  for  TGF  testing. 


Location 

Feed  Tube  ID  (inches) 

Shroud  Tube  ID  (inches) 

Stagnation  Chamber 

0.068 

0.40 

Stagnation  Chamber 

0.055 

0.33 

Test  Section 

0.068 

0.17 

Figure  4.27:  Cross-plane  PIV  setup  at  the  TGF  with  approximate  fields  of  view  for 
each  camera. 

At  the  conditions  in  Table  3.2,  each  injector  was  tested  individually  to  observe 
particles.  The  goal  was  to  see  particles  injected  from  the  stilling  chamber  because  for 
actual  testing,  there  would  be  little  flow  disturbances  from  the  shroud  tube  seen  in 
the  test  section  whereas  the  injector  just  upstream  of  the  test  section  would  cause  a 
momentum  deficit  in  its  wake,  possibly  disturbing  the  flow  over  the  model  being  studied. 
If  no  particles  were  seen,  however,  this  injector  near  the  test  section  windows  would  be 
easily  observed  by  eye  to  check  for  problems.  Particles  were  observed  from  two  of  the 


57 


three  injectors,  though  the  larger,  curved  shroud  tube  provided  very  promising  results. 
This  injector  used  a  shorter,  3  inch  feed  tube  with  a  0.068  inch  ID.  Figure  4.28  is  a 
typical  image  of  the  right  side  of  the  test  section  looking  upstream  at  M  =  0.8  and  P 
=  250  psf  using  this  injector.  The  third,  smaller  stagnation  chamber  injector  did  not 
produce  particles  that  persisted  through  the  test  section.  Similar  results  were  observed 
at  all  speed  and  pressure  test  points. 
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Figure  4.28:  Image  of  seeding  from  curved  0.4  inch  ID  shroud  at  M=0.8,  P=250  psf. 


The  second  tube  injecting  into  the  stilling  chamber  was  a  simple  tube  with  a  0.33 
inch  ID  and  0.055  inch  ID  feed  tube,  identical  to  that  used  in  the  AFIT  tunnel.  Few, 
if  any  particles  were  observed  in  the  test  section  with  this  injection  configuration.  A 
typical  image  of  the  particles  from  this  injector  is  seen  in  Figure  4.29.  This  image  was 
taken  with  the  tunnel  at  M  =  0.8  and  P  =  250  psf,  which  was  the  highest  speed  and 
lowest  pressure  settings  used.  A  higher  speed  decreases  the  time  between  injection  and 
passage  through  the  test  section  for  individual  particles,  and  a  lower  pressure  increases 
the  expansion  ratio  from  the  liquid  CO2,  thus  this  condition  is  the  most  favorable  for 
particle  persistence  through  the  test  section.  Since  few  particles  were  visible  from  this 
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injector  at  this  tunnel  condition,  data  was  not  collected  with  this  injector  configuration 
for  the  remaining  test  conditions. 
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Figure  4.29:  Image  of  seeding  from  simple  0.33  inch  ID  shroud  at  M=0.8,  P=250  psf 
showing  lack  of  particles. 

The  third  injector  was  installed  13  inches  upstream  of  the  laser  sheet  and  projected 
4  inches  into  the  tunnel  from  the  test  section  wall.  It  used  a  0.17  inch  ID  shroud  tube  and 
a  0.055  inch  feed  tube.  Because  of  its  close  proximity  to  the  laser  sheet  and  the  viewing 
windows,  the  particles  produced  were  easily  observed  by  both  the  PIV  system  and  by 
eye.  However,  the  close  proximity  also  meant  that  the  particles  did  not  spread  evenly 
across  the  FOV  of  either  camera,  so  discrete  particles  were  hard  to  discern  because  of 
the  large  concentration  over  a  very  small  region.  Figure  4.30  is  a  typical  image  of  the 
particles  produced  by  the  test  section  injector  and  Figure  4.31  overlays  this  image  on  the 
approximate  field  of  view  of  the  camera.  The  bloom  is  indicative  of  a  jet  in  a  crossflow 
profile  where  the  larger  particles  are  ejected  further  into  the  stream  and  smaller  particles 
respond  more  quickly  to  the  flow  change.  Since  the  main  goal  was  to  observe  particles 
injected  into  the  stilling  chamber  and  the  curved  injector  proved  effective  in  achieving 
this,  the  test  section  injector  and  the  other  stilling  chamber  injector  were  not  used  for 
the  remaining  test  conditions. 

4-3.2  Streamwise  Plane.  With  successful  particle  generation  using  the  large 
ID  feed  tube  and  the  large  ID,  curved  shroud  tube,  the  TGF  was  reconfigured  with  two 
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Figure  4.30:  Image  of  seeding  from  simple  0.17  inch  ID  shroud  at  M=0.8,  P=250  psf 
showing  typical  jet  in  a  crossflow  profile. 


Figure  4.31:  Approximate  PIV  image  location  within  the  test  section. 
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more  injectors  identical  to  the  first,  all  injecting  into  the  stilling  chamber.  The  PIV 
system  was  reconfigured  such  that  the  laser  sheet  was  parallel  with  the  flow  direction, 
but  at  an  angle  from  the  tunnel  walls.  This  was  necessary  to  allow  a  perpendicular  angle 
between  the  camera  and  laser  sheet.  This  setup  is  shown  in  Figure  4.32.  The  tunnel 
was  run  at  the  same  conditions  as  shown  in  Table  3.2  with  all  three  injectors  flowing 
at  each  test  point.  One  hundred  image  pairs  were  collected  at  a  frequency  of  10  Hz  for 
each  condition.  Figure  4.33  is  a  sample  image  collected  in  this  configuration.  Several 
reflections  are  present  in  the  left  side  of  the  image,  thus  this  data  was  removed  from  the 
figures  showing  the  results  for  all  the  conditions.  The  red  box  in  the  image  indicates 
the  area  that  PIV  figures  are  drawn  from.  Also  visible  in  this  image  is  the  variation  in 
particle  size  seen  in  the  FOV.  Many  large  particles  were  present  in  the  images.  Figure 
4.34  projects  the  PIV  results  presented  on  the  approximate  location  in  the  wind  tunnel 
test  section.  To  align  with  standard  convention,  the  data  has  been  rotated  so  that  the 
flow  is  from  left  to  right  in  each  vector  map. 


Figure  4.32:  Setup  for  PIV  data  collection  in  the  TGF. 


Figures  4.35  and  4.36  are  comparisons  of  the  three  different  test  condition  speeds 
at  approximately  the  same  stagnation  pressure.  In  these  and  the  remaining  figures, 
the  images  have  been  rotated  so  that  the  flow  direction  is  left  to  right  as  indicated  by 
the  vector  maps.  Additionally,  the  full  frame  of  the  camera  FOV  has  been  cropped  to 
remove  erroneous  data  that  was  caused  by  strong  reflections  in  the  images.  The  TGF 
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Figure  4.33:  Sample  PIV  image  collected  in  the  TGF.  The  red  box  indicates  area 

processed  for  PIV. 


Figure  4.34:  Sample  PIV  vector  map  overlayed  on  the  approximate  field  of  view  for 
the  PIV  arrangement  a  the  TGF. 
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operated  at  a  stagnation  temperature  between  293  K  and  300  K  so  the  actual  velocity 
using  equation  2.8  can  be  calculated.  Based  on  tunnel  data  collected,  the  approximate 
velocities  at  Mach  0.3,  0.5  and  0.8  were  100  m/s,  170  m/s  and  260  m/s,  respectively. 

Figure  4.37  is  the  streamwise  velocity  component  at  Mach  0.89  and  255  psf.  Figure 
4.38  is  at  Mach  0.79  and  501  psf.  At  the  lower  pressure  there  were  larger  variations  in 
speed  than  at  the  higher  pressure  condition  where  the  held  is  much  more  uniform.  This 
was  likely  due  to  a  variation  in  tunnel  conditions  during  data  collection  that  was  not 
known  at  the  time.  It  was  noted  following  data  collection  that  the  TGF  conditions  were 
not  stable  and  the  Mach  number  varied,  which  could  account  for  the  speed  variation  in 
the  data.  Each  vector  was  the  average  of  only  about  10  -  20  %  of  the  image  pairs  so 
some  vectors  could  be  the  result  from  the  beginning  of  the  data  collection  whereas  its 
neighboring  vectors  may  have  been  the  result  of  data  from  the  end  of  the  data  collection. 
For  the  same  reason,  the  turbulence,  as  indicated  by  the  RMS  fluctuations  normalized 
by  the  streamwise  velocity,  was  higher  than  expected,  as  seen  in  Figures  4.39,  4.40,  4.41 
and  4.42. 
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Figure  4.35:  Streamwise  velocity  maps  at  Pq  =  500  psf  and  varying  speeds. 
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Figure  4.36:  Vertical  velocity  maps  at  Pq  =  500  psf  and  varying  speeds. 


65 


St  ream  wise  Distance  (mm) 


Figure  4.37: 


Streamwise  velocity  profile  at  M  =  0.89,  P0  =  255  psf,  and  T0  =  294  K. 


Figure  4.38:  Streamwise  velocity  profile  at  M  =  0.79,  P0  =  501  psf,  and  T0  =  296  K. 
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Figure  4.39:  Streamwise  turbulence  profile  at  M  =  0.89,  P0  =  255  psf,  and  T0  =  294 
K  normalized  by  the  calculated  freestream  velocity,  286  m/s. 
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Figure  4.40:  Vertical  velocity  turbulence  profile  at  M  =  0.89,  Po  =  255  psf,  and  T0  = 
294  K  normalized  by  the  calculated  freestream  velocity,  286  m/s. 
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Figure  4.41:  Streamwise  velocity  turbulence  profile  at  M  =  0.79,  Po  =  501  psf,  and 
T0  =  296  K  normalized  by  the  calculated  freestream  velocity,  257  m/s. 
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Figure  4.42:  Vertical  velocity  turbulence  profile  at  M  =  0.79,  Po  =  501  psf,  and  T0  = 
296  K  normalized  by  the  calculated  freestream  velocity,  257  m/s. 
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The  turbulence  figures,  however,  are  the  result  of  only  100  image  pairs.  Therefore, 
the  turbulence  results  are  presented  as  a  rough  estimate  rather  than  a  statistically  valid 
result.  If  each  valid  vector  for  a  test  condition  is  assumed  to  be  an  independent  measure 
of  the  freestrearn  velocity,  however,  then  they  can  be  combined  into  a  single  ensemble 
of  several  thousand  data  points  versus  100.  Table  4.2  presents  these  results,  providing 
the  number  of  vectors  used  at  each  test  condition,  the  average  horizontal  and  vertical 
velocity  components,  and  the  normalized  fluctuations.  The  Mach  0.3,  P  =516  psf  and 
the  Mach  0.8,  P=255  psf  cases  have  the  largest  fluctuations,  but  these  are  explained 
above  for  the  high  speed  case  and  in  later  paragraphs  for  the  lower  speed  case. 


Table  4.2:  Turbulence  estimates  based  on  all  valid  vectors  at  each  test  condition. 


Mach 

P  (psf) 

#  Vectors 

U  avg 

|vu 

U’/U 

V’/U 

0.30 

517 

2931 

94.0 

4.5 

0.0186 

0.0011 

0.30 

1007 

14877 

107.4 

3.3 

0.0080 

0.0003 

0.30 

1506 

14955 

109.9 

3.2 

0.0080 

0.0003 

0.5 

500 

12630 

181.7 

7.4 

0.0086 

0.0004 

0.5 

1000 

15044 

175.3 

3.7 

0.0080 

0.0002 

0.5 

1500 

15005 

173.4 

3.0 

0.0080 

0.0002 

0.89 

255 

2282 

272.7 

11.1 

0.0203 

0.0009 

0.79 

501 

14107 

258.6 

8.6 

0.0082 

0.0003 

Figures  4.43,  4.44  and  4.45  are  the  streamwise  velocity  maps  for  a  nominal  Mach 
0.5  velocity  and  pressures  of  500,  1000  and  1500  psf.  The  actual  tunnel  data  for  these 
three  events  was  not  available,  so  true  velocity  comparisons  could  not  be  made.  The 
calculated  velocity  at  Mach  0.5  is  approximately  170  m/s  but  at  the  lower  pressures, 
PIV  data  indicates  an  approximately  8%  higher  velocity  than  this.  However,  if  the 
actual  speed  was  higher  during  the  lower  pressure  runs  because  the  TGF  was  set  to  a 
higher  velocity,  then  the  PIV  data  may  be  accurate.  Vertical  velocity  remains  low  and 
freestrearn  turbulence  is  less  than  2%  for  each  case,  as  shown  in  Figures  4.46,  4.48,  4.50, 
4.47,  4.49  and  4.51. 
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4.43: 


Streamwise  velocity  profile  at  M  =  0.5,  P0  =  500  psf,  and  T0  =  294  K. 
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Figure  4.44: 


Streamwise  velocity  profile  at  M  =  0.5,  P0  =  1000  psf,  and  T0  =  294  K. 
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Figure  4.45: 


Streamwise  velocity  profile  at  M  =  0.5,  P0  =  1500  psf,  and  T0  =  294  K. 
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Figure  4.46:  Streamwise  turbulence  profile  at  M  =  0.5,  Po  =  500  psf,  and  To  =  294  K 
normalized  by  the  calculated  freestream  velocity,  170  m/s. 


71 


-10 

S'20 

£-30 

> 

-40 

-20  0  20  40  60 

Streamwise  Direction  (mm) 


I  I  I  I  I  I  I  I  I  I  I  I  I  I  L 


Vy_prime_normalized 


0.02 

0.015 

0.01 

0.005 

0 


Figure  4.47:  Vertical  velocity  turbulence  profile  at  M  =  0.5,  Po  =  500  psf,  and  T0  = 
294  K  normalized  by  the  calculated  freestream  velocity,  170  m/s. 
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Figure  4.48:  Streamwise  turbulence  profile  at  M  =  0.5,  P0  =  1000  psf,  and  T0  =  294 
K  normalized  by  the  calculated  freestream  velocity,  170  m/s. 
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Figure  4.49:  Vertical  velocity  turbulence  profile  at  M  =  0.5,  P0  =  1000  psf,  and  T0  = 
294  K  normalized  by  the  calculated  freestream  velocity,  170  m/s. 
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Figure  4.50:  Streamwise  turbulence  profile  at  M  =  0.5,  P0  =  1500  psf,  and  T0  =  294 
K  normalized  by  the  calculated  freestream  velocity,  170  m/s. 
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Figure  4.51:  Vertical  velocity  turbulence  profile  at  M  =  0.5,  Po  =  1500  psf,  and  T0  = 
294  K  normalized  by  the  calculated  freestream  velocity,  170  m/s. 
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Figures  4.52,  4.53  and  4.54  are  the  streamwise  velocity  maps  for  a  nominal  Mach 
0.30  velocity  and  pressures  of  500,  1000  and  1500  psf.  Based  on  tunnel  instrumentation, 
the  expected  velocities  for  these  three  cases  were  approximately  100  m/s,  104  m/s  and 
104  m/s.  The  graphs  are  centered  on  these  expected  values  and  color  contours  are 
placed  for  the  range  ±  15  m/s.  Following  data  collection  at  the  500  psf  condition  we 
found  that  at  least  one  CO2  injector  had  clogged,  so  the  mass  flow  had  been  reduced, 
decreasing  the  number  of  particles  and  decreasing  the  smoothness  of  that  data.  The  PIV 
data  for  the  three  conditions  showed  average  streamwise  velocities  that  closely  matched 
expected  velocities.  Streamwise  turbulence  remained  below  1.5%  across  much  of  the 
field  whereas  the  vertical  velocity  component  turbulence  was  less  than  1%  in  the  the  two 
higher  pressure  cases.  The  low  pressure  case  had  pockets  of  turbulence  near  2%,  but 
again  this  is  the  test  case  that  experienced  a  reduction  of  particles,  so  a  decrease  in  the 
smoothness  is  expected.  Turbulence  profiles  are  given  in  Figures  4.55,  4.57,  4.59,  4.56, 
4.58  and  4.60. 
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Figure  4.52:  Streamwise  velocity  profile  at  M  =  0.30,  P0  =  517  psf,  and  T0  =  298  K. 


4-3.3  Sources  of  Error  in  AFRL  Wind  Tunnel  Testing.  Two  sources  of  error 
have  been  discussed  above,  causing  less  uniform  results  in  the  PIV  data:  first,  the  fluc¬ 
tuation  in  tunnel  conditions  during  data  collection,  and  second,  the  clogging  of  the  C02 
lines.  In  addition  to  these,  large  particles  present  in  the  flow  and  strong  reflections  both 
caused  uncertainty  in  the  data.  The  large  particles  caused  pixel  saturation  with  streak- 
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Figure  4.53:  Streamwise  velocity  profile  at  M  =  0.30,  P0  =  1007  psf,  and  T0  =  298  K. 


Figure  4.54:  Streamwise  velocity  profile  at  M  =  0.30,  P0  =  1506  psf,  and  T0  =  300  K. 
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Figure  4.55:  Streamwise  turbulence  profile  at  M  =  0.30,  Po  =  517  psf,  and  T0  =  298 
K  normalized  by  the  calculated  freestream  velocity,  103  m/s. 


76 


Vy_p  rime_n  orm  alize  d 

0.02 

0.015 

0.01 

0.005 

0 

Figure  4.56:  Vertical  velocity  turbulence  profile  at  M  =  0.30,  P0  =  517  psf,  and  T0  = 
298  K  normalized  by  the  calculated  freestrea.ru  velocity,  103  m/s. 
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Figure  4.57:  Streamwise  turbulence  profile  at  M  =  0.30,  P0  =  1007  psf,  and  T0  =  298 
K  normalized  by  the  calculated  freestream  velocity,  104  m/s. 
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Figure  4.58:  Vertical  velocity  turbulence  profile  at  M  =  0.30,  P0  =  1007  psf,  and  T0 
=  298  K  normalized  by  the  calculated  freestream  velocity,  104  m/s. 
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Figure  4.59:  Streamwise  turbulence  profile  at  M  =  0.30,  Po  =  1506  psf,  and  T0  =  300 
K  normalized  by  the  calculated  freestream  velocity,  104  m/s. 
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Figure  4.60:  Vertical  velocity  turbulence  profile  at  M  =  0.30,  Po  =  1506  psf,  and  T0 
=  300  K  normalized  by  the  calculated  freestream  velocity,  104  m/s. 

ing,  decreasing  the  number  of  valid  vectors  available  to  render  the  figures  above.  Also, 
the  reflections  present  in  the  image  had  the  same  effect,  reducing  the  area  used  for  PIV 
to  about  a  third  of  the  image  field  of  view.  Finally,  as  mentioned  above,  these  results 
are  based  on  a  single  set  of  100  image  pairs  for  each  condition.  Though  particle  coverage 
was  sufficient  to  compute  these  results,  typically  only  about  10-20  %  of  available  vectors 
were  used  to  determine  the  mean  vector  magnitudes  presented  in  these  results.  Thus, 
more  particles  and  more  image  pairs  could  improve  the  results. 
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V.  Conclusions 


5.1  Overview  of  Research  Effort 

PIV  is  a  highly  effective  measurement  technique,  but  its  usefulness  depends  on 
careful  selection  of  particles  for  flow  seeding.  Particle  characteristics  such  as  size  and 
density  are  important  in  determining  the  response  of  the  particle  to  the  flow  dynamics. 
Additionally,  particle  size  affects  the  signal  recorded  by  cameras  used  to  collect  PIV  data. 
In  many  circumstances,  particularly  in  closed-circuit  tunnels,  the  particles  can  affect  the 
safety  and  performance  of  the  facility  where  they  are  to  be  used.  Facilities  such  as  the 
AFRL/RB  TGF  have  difficulty  using  PIV  because  of  the  impact  of  the  seed  material  on 
tunnel  components.  Seeding  materials  can  accumulate  within  the  tunnel,  either  attaching 
to  wall  surfaces  or  system  components  necessitating  additional  maintenance  for  tunnel 
operation.  Thus,  a  method  of  cleanly  seeding  the  tunnel  for  PIV  measurements  is  sought. 

To  this  end,  research  has  been  performed  on  the  introduction  of  solid  CO2  particles 
for  the  application  of  PIV  seeding.  These  dry  ice  particles  are  generated  by  allowing  liquid 
CO2  to  expand  through  an  orifice  into  a  shroud  tube  which  leads  to  solidified  particles 
which  are  injected  into  the  wind  tunnel.  The  present  study  comprised  determining  the 
effects  of  changing  the  orifice  tube  and  shroud  tube  size  on  particle  size  by  using  a  Malvern 
Spraytec  particle  size  analyzer  to  obtain  mean  particle  sizes.  In  addition  to  the  varying 
sizes  of  tubes,  a  different  shroud  tube  using  a  Kenics®  static  mixer  and  heat  exchanger 
tube  was  studied.  Two  combinations  of  injection  methods  were  assessed  in  the  AFIT  6 
inch  by  6.5  inch  supersonic  blowdown  wind  tunnel  at  Mach  1.92.  A  10°  ramp  in  the  test 
section  led  to  an  oblique  shock  wave,  and  particle  response  was  characterized.  In  this 
study  three  separate  injection  ports  with  identical  orifice  and  shroud  tube  arrangements 
were  used  in  this  tunnel.  This  significantly  increased  particle  population  in  the  test 
section,  which  generally  improves  the  spatial  resolution  and  accuracy  of  the  PIV  data. 
Also  three  injection  methods  were  qualitatively  observed  in  the  AFRL  TGF  at  varying 
subsonic  conditions  to  demonstrate  the  scaling  effects  and  injection  geometry  of  a  larger 
tunnel  (24  inch  by  24  inch  test  section).  Finally,  PIV  data  was  collected  using  the 
preferred  injection  method  in  the  TGF  at  varying  subsonic  conditions,  demonstrating 
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the  applicability  of  the  clean  seeding  method  for  PIV  measurements  at  Mach  numbers 
of  0.3,  0.5  and  0.8. 

5. 2  Conclusions 

Particle  size  data  collected  using  the  Malvern  Spraytec  suggests  that  the  static 
mixing  shroud  tube  produced  particles  with  a  Sauter  mean  diameter,  D32,  of  approxi¬ 
mately  3  times  smaller  than  that  for  a  simple  straight  shroud  tube  of  nearly  equal  size. 
This  tube  was  also  more  likely  to  clog,  causing  a  reduction  of  particles.  In  separate 
experiments,  the  two  different  injector  shroud  types  were  then  used  in  the  AFIT  6  inch 
by  6.5  inch  wind  tunnel  and  particle  response  was  assessed  across  an  oblique  shock  wave. 
Three  injectors  of  each  type  were  employed,  injecting  into  the  stagnation  chamber  of  the 
tunnel.  The  expectation  was  that  smaller  particles  would  be  produced  and  they  would 
respond  to  the  oblique  shock  in  a  shorter  distance  than  larger  particles.  The  raw  images 
suggested  that  smaller  particles  were  being  produced  when  the  shroud  with  the  mixer 
was  used  as  it  was  necessary  to  increase  the  laser  power  and  open  the  lens  aperture  to 
obtain  images  of  similar  quality  to  those  captured  using  the  simple  shroud  tube.  This 
shroud  tube  increased  the  likelihood  of  clogging,  however,  causing  a  reduction  in  particles 
when  clogging  occurred. 

The  particles  generated  using  the  simple  0.33  inch  shroud  tube  and  0.055  inch  feed 
tube  closely  matched  the  freestream  horizontal  velocity  expected  from  tunnel  instrumen¬ 
tation  within  3%  using  the  larger,  64  x  64  interrogation  region,  which  corresponded  to 
2.43  mm  by  2.43  mm.  The  horizontal  velocity  far  downstream  of  the  shock  was  within 
approximately  3%  of  the  calculated  values  as  well.  The  vertical  velocity  component  of 
the  PIV  results  was  15  m/s  upstream  and  65  m/s  downstream  of  the  shock.  The  expected 
vertical  velocity  component  upstream  of  the  shock  was  assumed  to  be  zero,  and  based 
on  tunnel  instrumentation  and  isentropic  relations,  the  component  downstream  of  the 
shock  was  calculate  to  be  77  m/s.  The  discrepancy  in  the  vertical  component  suggests 
about  1°  -  2°  of  flow  angularity. 

Using  the  static  mixing  shroud  tubes  and  same  feed  tubes  for  similar  test  condi¬ 
tions,  the  differences  between  the  computed  and  experimental  velocities  upstream  and 
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downstream  of  the  shock  were  also  low.  There  was  virtually  no  effect  on  the  shroud  tube 
on  these  mean  velocity  values.  In  both  cases,  using  the  smaller  interrogation  region  (IR) 
caused  a  bias  towards  lower  velocities  which  is  in  accord  with  the  literature  and  the  soft¬ 
ware  manufacturer  specifications.  A  reason  for  exploring  a  smaller  IR  was  to  determine 
whether  the  measured  response  of  the  particles  to  the  shock  could  be  improved  in  this 
experiment.  There  was  no  improvement  when  going  to  the  smaller  IR  in  either  data  set, 
suggesting  that  the  IR  size  was  not  a  factor  in  the  response  calculations. 

Across  the  shock,  the  particles  recovered  to  expected  velocities  based  on  tunnel 
instrumentation  within  approximately  15  mm  normal  to  the  shock.  This  would  suggest 
particle  sizes  of  approximately  2  /mi  for  both  data  sets,  following  the  methodology  laid 
out  by  Melling  [13].  Schlieren  images  showed  a  strong  interaction  between  the  oblique 
shock  and  the  tunnel  floor  boundary  layer  effects;  shock  motion  of  approximately  4  mm 
was  observed  in  the  freestream.  Particle  response  was  calculated  based  on  the  PIV  data 
which  is  an  array  of  mean  velocity  data.  Therefore  transient  shock  conditions  would  be 
averaged  during  the  data  processing,  causing  a  smeared  shock  region  in  the  velocity  data 
results.  Thus,  the  assumption  that  the  shock  acted  as  a  step  function  is  not  necessarily 
valid  for  this  setup.  This  shock  motion  suggests  the  response  comparison  to  a  theoretical 
2  /mi  particle  is  likely  conservative  and  actual  particle  response  could  more  closely  match 
that  of  smaller  particles  using  Mclling’s  method. 

For  the  first  time,  the  C02  seeding  method  was  scaled  up  and  applied  to  the 
TGF  successfully.  By  positioning  the  laser  sheet  perpendicular  to  the  flow  direction, 
observations  could  be  made  on  whether  particles  injected  into  the  stagnation  region 
were  persisting  long  enough  to  pass  through  the  test  section  in  this  tunnel.  The  0.068 
inch  ID  feed  tube  and  0.40  inch  ID  curved  shroud  to  direct  particles  downstream  from 
the  stagnation  chamber  produced  the  best  particle  coverage  so  this  combination  was 
replicated  two  times  to  provide  three  identical  injectors  at  three  points  on  either  side 
of  the  stagnation  chamber.  A  simple  shroud  tube  injecting  in  the  stagnation  region 
normal  to  the  streamwise  direction  was  ineffective.  The  curved  shroud  tube  was  critical 
to  the  success  because  the  particles  were  directed  toward  the  test  section,  reducing  the 
initial  velocity  difference  between  the  particles  and  the  freestream.  After  repositioning 
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the  laser  sheet  and  camera,  PIV  measurements  were  taken  at  three  subsonic  speeds  at 
varying  pressures.  In  all  cases,  the  seeding  was  sufficient  to  provide  mean  velocities  in  the 
flow  field  which  matched  predicted  values  within  a  few  percent.  Freestream  turbulence 
was  typically  between  1-2%  except  where  problems  with  particle  coverage  due  to  clogging 
occurred  and  when  tunnel  settings  drifted  during  data  collection.  The  combination  of 
injectors  used  in  the  TGF  produced  some  large  particles  that  saturated  camera  pixels 
introducing  error  into  the  IR  correlations.  Collecting  and  analyzing  100  image  pairs 
using  a  range  filter  limited  the  impact  of  the  large  particles  when  computing  mean 
velocities,  though  improved  methods  of  dealing  with  the  large  particles  are  needed  so 
that  the  method  may  be  employed  robustly.  Importantly,  no  residue  persisted  in  the 
tunnel  following  the  data  collection  using  the  C02  to  seed  the  flow. 

5.3  Impact  of  Research 

For  the  first  time  C02  was  used  to  seed  the  TGF  wind  tunnel  for  PIV.  Facilities 
such  as  the  TGF  can  benefit  from  this  method  because  of  the  low  cost  and  availability 
of  liquid  C02,  the  ability  to  influence  particle  size  by  injector  design,  and  because  the 
seeding  does  not  cause  maintainability  or  safety  problems.  The  characterization  of  the 
particle  response  may  lead  to  improvements  by  providing  a  means  of  calculating  the 
particle  lag  in  the  fluid.  This  is  necessary  to  draw  meaningful  velocity  data  from  PIV 
in  dynamic  fluid  motion  and  to  resolve  small-scale  motion  in  areas  of  interest  where  the 
relaxation  distance  of  the  particle  is  much  closer  in  magnitude  to  the  scale  being  resolved. 
Finally,  the  effect  on  particle  size  due  to  varying  the  injection  method  may  be  used  to 
size  the  particles  for  the  geometry  of  the  tunnel  it  is  used  in.  The  distance  between  the 
injection  point  and  the  test  section  as  well  as  the  speed  of  the  air  flow  of  the  tunnel  are 
important  for  determining  the  residence  time  the  particles  must  be  able  to  achieve  to  be 
useful  before  sublimating  completely. 

5.4  Future  Work 

This  and  previous  research  have  demonstrated  the  capabilities  and  limitations  of 
various  injection  methods  and  injection  locations,  and  have  successfully  scaled  the  in- 
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jection  method  for  larger  wind  tunnel  use.  Future  work  may  be  devoted  to  acquiring 
particle  response  data  for  a  stationary  shock  that  is  not  influenced  by  boundary  layer 
interaction.  By  designing  a  vertical  wedge  rather  than  a  ramp,  the  shock-boundary  layer 
interaction  can  be  minimized  and  a  shock  attached  to  the  wedge  can  be  observed  in  the 
horizontal  plane.  Also,  a  method  for  limiting  particle  dispersion  to  enhance  the  concen¬ 
tration  of  particles  in  a  desired  area  would  enhance  the  accuracy  of  mean  velocities  and 
enable  smaller  interrogation  regions,  as  well  as  make  time-resolved  PIV  possible.  In  this 
and  previous  research,  the  particles  have  mainly  been  injected  perpendicular  to  the  flow 
direction  within  the  stagnation  region.  This  method  has  been  suitable  for  observing  an 
oblique  shock  wave  as  well  as  boundary  layer  flow  in  previous  work  because  it  provides  an 
even  distribution  of  particles  across  the  entire  test  section  cross-sectional  area.  However, 
in  larger  tunnels,  an  even  distribution  can  be  unnecessary.  Given  that  PIV  often  focuses 
on  a  small  subregion  of  a  large  tunnel,  the  ability  to  inject  particles  that  cover  just  this 
region  could  enhance  the  ability  to  conduct  PIV  without  having  to  utilize  an  excessive 
amount  of  C02  to  seed  the  flow.  In  larger  facilities  such  as  the  TGF,  a  structure  could 
be  installed  completely  within  the  stagnation  region  such  that  the  injector  tubes  could 
be  aligned  with  the  flow  direction  as  well  as  translated  vertically  and  horizontally  such 
that  the  particle  concentration  could  be  localized  rather  than  dispersed.  These  areas  of 
research  may  be  unique  to  the  tunnel  using  this  technique,  so  future  work  could  include 
the  full  scale  introduction  and  use  of  dry  ice  seeding  in  a  facility  such  as  the  TGF  for  a 
specific  test. 
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Appendix  A.  Theoretical  Response  Curves  Using  Melling’s  Approach 

The  following  spreadsheet  is  used  to  calculate  the  theoretical  response  of  particles  as 
proposed  by  Melling.  The  variables  are  calculated  using  the  following  expressions,  given 
an  upstream  stagnation  temperature,  T0j,  and  Mach  number,  Mi,  as  well  as  the  shock 
angle,  /3.  In  all  cases,  the  subscript,  “1”  refers  to  upstream  of  the  shock,  and  “2”  refers 
to  downstream  of  the  shock. 


Ui  =  Mu 


7-R 


Toa 


\  + 


Where  7  =  1.4,  and  R  =  287 JKg  lK  1. 

To  ensure  the  f3  used  results  with  the  correct  ramp  angle,  9,  can  be  calculated: 


tan  9  =  2  cot  /3 


M(  sin  (3  -  1 
Ml  (7  +  cos  2/3)  +  2 


The  Mach  number  normal  to  the  shock  is: 


=  Mi  sin  /3 


The  normal  velocity  upstream  of  the  shock  is: 


VnA  =  Ui  sin/3 


The  total  temperature  downstream  of  the  shock  is: 


T2  =  7i  [2  +  (7  -  1)  Ml,} 


27  Via  -  (7  ^  1) 
(7  +  1  ?Ml, 


The  Mach  number  normal  to  and  downstream  of  the  shock  is: 


2  1  +(2=i)Af21 

Ml 2  =  V  2  ;  n,i 


7 Mil  -  ¥ 
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The  velocity  normal  to  and  downstream  of  the  shock  is: 


14,2  —  2  \/  1RT2 


The  Mach  number  downstream  of  the  shock  is: 


M2 


d4i,2 

sin  (/3  -  9) 


The  velocity  downstream  of  the  shock  is: 


u2  =  m2S/1rt  cos  9 


Finally,  Mclling’s  relations  can  be  applied  to  determine  the  theoretical  particle  response: 


=  18/x 
-  Pv% 


Calculations  for  Particle 


T1  = 

265 

Ml  = 

1.92 

U1  = 

475  3289018 

Beta  = 

41 

Theta  = 

0  175890805 

Mnl  = 

1  259633336 

Vnl  = 

311.8438178 

T2  = 

308  8390331 

Mn2  = 

0  651711701 

Vn2  = 

229  5760486 

M2  = 

1  2682333 

U2  = 

439  8628391 

10  07780077  deg 


Response 

p=1180  kg/m3 


dp 

c 

0  0000002 

2  82263E-07 

0  0000004 

1  12905E-06 

0  0000006 

2.54037E-06 

0  0000008 

4  51621E-06 

0  000001 

7  05657E-06 

0  0000012 

1  01615E-05 

0  0000014 

1  38309E-05 

0  0000016 

1  80648E-05 

0  0000018 

2  28633E-05 

0  000002 

2  82263E-05 

0  0000022 

3.41538E-05 

0  0000024 

4  06459E-05 

0  0000026 

4  77024E-05 

0  0000028 

5.53235E-05 

0  000003 

6  35091E-05 
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Theoretical  Particle  Responses  Using  AFIT  Wind  Tunnel  Conditions 


Time(s) 

Vn  I.Oum 

Dist.  Traveled  (mm) 

Vn  2.0um 

Dist  Traveled 

Vn  3.0um 

Distance  (mm) 

311  8438178 

-0  602665882 

311  8438178 

-0  618010362 

311.8438178 

-0.621127227 

311  8438178 

-0.30640885 

311.8438178 

-0.31041204 

311.8438178 

-0  311201206 

311  8438178 

0 

311.8438178 

0 

311.8438178 

0 

0 

311  8438178 

0.30640885 

311.8438178 

0  31041204 

311.8438178 

0  311201206 

0.000001 

300.9738832 

0.602665882 

308.9802625 

0.618010362 

310.558594 

0.621127227 

0  000002 

291  540179 

0.89011244 

306.2163812 

0  922892905 

309.2934485 

0  929797985 

0  000003 

283  3529381 

1.169912643 

303  5487043 

1.225154198 

308  0480677 

1  237233091 

0.000004 

276.2474668 

1.443076792 

300.9738832 

1.524885483 

306.8221428 

1.543451847 

0  000005 

270.0808323 

1.710481702 

298  4886858 

1.822174822 

305.6153699 

1  848473257 

0  000006 

264  7289873 

1.972888333 

296  0899925 

2  117107214 

304.4274497 

2  152316025 

0.000007 

260.0842746 

2.230957101 

293.7747923 

2.4097647 

303.2580877 

2.454998566 

0.000008 

256.0532619 

2  485261163 

291.540179 

2.700226463 

302.106994 

2.756539005 

0  000009 

252.5548618 

2  736297944 

289.3833476 

2  988568933 

300  9738832 

3  056955183 

0.00001 

249.5187007 

2.984499146 

287  3015906 

3.274865875 

299.8584743 

3.356264666 

0.000011 

246.8837037 

3.230239431 

285.2922948 

3  559188492 

298.7604909 

3.654484742 

0  000012 

244  5968655 

3.473843956 

283  3529381 

3  841605504 

297.6796607 

3  95163243 

0.000013 

242.6121842 

3.715594916 

281.481086 

4.122183241 

296.6157157 

4.247724484 

0.000014 

240.8897364 

3.955737221 

279.6743889 

4.400985725 

295.5683921 

4  542777395 

0  000015 

239  3948733 

4  19448342 

277  9305788 

4  678074748 

294  5374302 

4  836807397 

0.000016 

238.0975246 

4.432017979 

276.2474668 

4  953509951 

293.5225745 

5  129830471 

0.000017 

236.971593 

4  66850099 

274.6229401 

5.227348901 

292.5235734 

5  421862348 

0  000018 

235  9944292 

4.904071393 

273  0549596 

5499647159 

291  540179 

5.712918511 

0  000019 

235.1463768 

5.13884977 

271.5415569 

5.770458354 

290.5721477 

6  003014204 

0.00002 

234.4103765 

5.372940769 

270.0808323 

6.039834246 

289.6192395 

6  292164433 

0  000021 

233.771623 

5.606435214 

268  6709522 

6.307824796 

288.6812179 

6  580383967 

0  000022 

233.2172671 

5.839411927 

267  3101469 

6  574478223 

287.7578506 

6  867687347 

0.000023 

232.7361576 

6.071939314 

265.9967082 

6.839841071 

286.8489085 

7.154088884 

0.000024 

232.3186164 

6.304076744 

264.7289873 

7.103958261 

285.9541663 

7.439602669 

0.000025 

231  9562445 

6.535875742 

263  5053928 

7  366873152 

285.0734022 

7.724242569 

0.000026 

231.6417522 

6.767381025 

262.324389 

7.628627593 

284.2063978 

8.008022237 

0.000027 

231.3688135 

6  998631401 

261.1844932 

7.889261977 

283.3529381 

8  290955111 

0  000028 

231  1319377 

7.22966055 

260  0842746 

8  14881529 

282  5128115 

8  573054422 

0.000029 

230.9263601 

7.460497702 

259.0223522 

8.407325163 

281.6858098 

8  854333191 

0.00003 

230.7479452 

7.691168227 

257.9973929 

8.664827914 

280.8717279 

9  134804237 

0.000031 

230.593104 

7  92169414 

257.00811 

8  9213586 

280  0703639 

9414480178 

0.000032 

230.4587218 

8  152094548 

256.0532619 

9.176951056 

279.2815191 

9.693373437 

0.000033 

230.3420953 

8  382386035 

255.1316499 

9  43163794 

278.5049981 

9.97149624 

0.000034 

230.2408785 

8  612582992 

254.2421172 

9.685450772 

277.7406082 

10  24886062 

0  000035 

230.1530354 

8  842697909 

253.383547 

9  938419976 

276.98816 

10.52547844 

0.000036 

230.0767989 

9.072741626 

252.5548618 

10  19057492 

276.2474668 

10.80136134 

0  000037 

230.0106354 

9  302723551 

251.7550212 

10  44194394 

275.5183451 

11  07652082 

0  000038 

229.953214 

9  532651848 

250.9830212 

10  6925544 

274  800614 

11  35096818 

0.000039 

229.9033796 

9.762533603 

250.2378928 

10.94243269 

274.0940957 

11.62471453 

0  00004 

229  8601298 

9.992374965 

249.5187007 

11  19160431 

273.3986149 

11  89777084 

0  000041 

229  8225945 

10  22218127 

248.824542 

11  44009386 

272  7139992 

12  17014788 

0.000042 

229.7900187 

10.45195715 

248.1545454 

11.68792507 

272.0400789 

12  44185626 

0  000043 

229  7617471 

10  68170663 

247.5078699 

11  93512085 

271.3766868 

12.71290643 

0  000044 

229  737211 

10  9114332 

246.8837037 

12.18170334 

270.7236586 

12  98330868 

0.000045 

229.7159168 

11.14113987 

246.2812633 

12.42769386 

270.0808323 

13.25307312 

0.000046 

229.6974362 

11.37082929 

245.6997925 

12.67311304 

269.4480485 

13  52220972 

0  000047 

229  6813974 

11  60050373 

245.1385614 

12  91798075 

268.8251503 

13  79072829 

0.000048 

229.6674778 

11  83016517 

244.5968655 

13  1623162 

268.2119833 

14.05863848 

0.000049 

229.6553974 

12  05981532 

244.0740248 

13  4061379 

267.6083955 

14  32594979 

0  00005 

229  6449131 

12  28945568 

243  5693829 

13  64946375 

267  0142372 

14  59267159 

0.000051 

229.6358141 

12.51908755 

243.0823066 

13  89231099 

266.4293611 

14  85881308 

5.2E-05 

229.6279174 

12.74871204 

242.6121842 

14  1346963 

265.8536222 

15  12438333 

5  3E-05 

229  621064 

12  97833013 

242  1584258 

14  37663574 

265.2868778 

15  38939127 

5  4E-05 

229.6151162 

13.20794267 

241.7204617 

14.61814484 

264.7289873 

15.65384567 

5.5E-05 

229.6099542 

13  43755038 

241.2977421 

14.85923858 

264.1798124 

15.91775518 

5  6E-05 

229  6054743 

13  66715391 

240  8897364 

15  09993142 

263.6392169 

16  18112832 

5  7E-05 

229  6015863 

13  89675381 

240.4959325 

15  3402373 

263.1070669 

16.44397347 

5.8E-05 

229.5982121 

14  12635056 

240.115836 

15  5801697 

262.5832304 

16.70629887 

5.9E-05 

229  5952836 

14.35594457 

239.7489698 

15  81974163 

262.0675774 

16  96811265 

6E-05 

229.5927421 

14  58553621 

239  3948733 

16  05896561 

261.5599803 

17  2294228 

6. IE-05 

229.5905364 

14  81512579 

239.0531022 

16.29785378 

261.060313 

17.49023718 

6.2E-05 

229.5886222 

15  04471358 

238.7232274 

16.53641781 

260.5684518 

17  75056355 

6  3E-05 

229  5869608 

15  27429982 

2384048348 

16  77466899 

260  0842746 

18  01040951 

6.4E-05 

229.585519 

15  50388471 

238.0975246 

17  01261821 

259.6076615 

18  26978259 

6.5E-05 

229.5842677 

15.73346844 

237.8009113 

17.25027597 

259.1384943 

18.52869017 

6.6E-05 

229  5831817 

15  96305115 

237  5146224 

1748765243 

258.6766565 

18  78713951 

6.7E-05 

229.5822392 

16  19263298 

237  2382986 

17.72475738 

258.2220339 

19  04513779 

6.8E-05 

229.5814212 

8.386864657 

236.971593 

9.66772322 

257.7745135 

10.28080433 
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